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Abstract 

PID controllers are widely used in industries due to 
simple functioning and robustness, yet the determination 
of controller gains is a big challenge more so when the 
plant to be control is multi input multi output (MIMO) 
and nonlinear in nature. The aim of this paper is to string 
together the qualities of conventional PID controller and 
a recent technique Genetic Algorithm (GA). In this work 
controller parameters are adjusted by means of fitness 
function that can evaluate the optimum PID parameters 
based on system error. Simulation study has been carried 
out for a turbo- generator and the results have been 
compared with decoupled control and diagonal recurrent 
neural network (DRNN) based PID tuning methods. A 
simulation result shows that the proposed control scheme 
results in smaller controller effort. 

Keywords: Genetic algorithm, PID Tuning, MIMO 
system 

1. Introduction 

PID controllers are widely used in industries due to their 
simple functioning and robustness but determination of 
controller gains has always been a big challenge. Most 
industrial processes are multivariable and their nonlinear 
and interactive behavior imposes difficulties in 
determination of the parameters of controller. In such a 
situation progressive technique like GA has been shown 
to be capable of high performance without experiencing 
the difficulties [1], [10]. GA offers an adaptive heuristic 
search algorithm based on the evolutionary ideas of 
selection and genetics [2]. It is inspired by Darwin’s 
theory about evolution- “Survival of fittest”. GA 
represents an intelligent exploitation of random search 
used to solve optimization problems. Takhiro and Sigeru 
[3] applied GA to tune the PID gains for a single input 
single output (SISO) water bath system and proved that 
real number encoding is simple to operate and better 
than bit string encoding method. Mitsukura, Yamamoto 
and Kaneda [4, 5] proposed a new tuning method of PID 
parameters based on the relationship between the PID 
control law and the generalized minimum variance 
control show the effectiveness of proposed method for 



the nonlinear SISO Hammerstein model. Maja and Kabra 
[8] presented a tuning method for multivariable controller 
that was based on decoupling control and GA. As this 
method tends to achieve the shortest direct paths through 
the system (without inverting the model), it usually 
results in low energy consumption which can further be 
directly controlled with pole placement. Jih-Gau Juang et 
al. [12] used GA to tune the PID parameters to control a 
laboratory helicopter called twin rotor system. All 
parameters of the controller are obtained by a real- value- 
type genetic algorithm (RGA) with an integral of time 
multiplied by the square error criterion. Jin-Sung Kim et. 
al [13] proposed tuning method which was based on 
improved GA. For better performance they defined new 
objective function in the sense of root mean square error. 
In addition they made improvement by replacing the 
uniform distribution random number generation in 
conventional GA to newly designed hybrid random 
generator composed of Cauchy distribution and linear 
congruential generator which provides independent and 
different random numbers at each individual steps in 
Genetic operations. In many research works GA has been 
widely applied to tune the PID parameters such as a 
nonlinear model of thermal process [9], gasifier problem 
[11], missile control system, DC motor position control, 
LTI system, liquid level system [14-17]. Monica 
Patrascu et al. [18] applied GA to tune the parameters of 
PID controller for nonlinear model MIMO system. They 
have focused on the initialization procedure of the 
algorithm that can influence the outcome of the 
optimization procedure and showed that roulette wheel 
selection is better choice for selection. Daniel Carmona 
Morales et.al [20] developed a Graphical user interface 
system for optimum tuning of multivariable PID 
controller. From the educational point of view, this tool 
provides students with the necessary means to 
consolidate their knowledge on these control structures. 
In this work GA is chosen to tune the parameters of PID 
controller for a MIMO Turbo- Generator plant which is 
highly nonlinear and has strong coupling among 
variables. The simulation results have obtained by the 
proposed approach been compared with those obtained 
by fixed decoupled PID configuration and also with 
DRNN based tuning of controller parameters [19]. 
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The rest of paper is organized as follows. In the next 
section basic steps of GA are described. In Section 3 
Genetic operator such as Selection, Crossover and 
mutation are illustrated. Section 4 describes the 
evaluation procedure based on fitness function and 
fitness value. In section 5 the basic ideas of the tuning 
method are presented. Section 6 describes the Turbo- 
Generator plant considered and simulation studies are 
presented in Section 7. Concluding remarks are in last 
Section. 

2. Genetic Algorithm 

GA is a stochastic search method that emulates the 
process of natural evolution. It starts with the minimal 
knowledge of the correct solution and depends entirely 
on responses from its environment and genetic operators 
e.g. selection or reproduction, crossover and mutation to 
get the optimal solution [1]. By starting at several 
independent points and searching in parallel, the GA 
avoids local minima and converges to sub optimal 
solutions. The GA is typically initialized with a random 
population consisting of between 20-100 individuals. 
This population is usually represented by a real valued 
number or a binary string called a chromosome. 
Individual performance is measured by the fitness 
function. The fitness function assigns each individual a 
corresponding number called its fitness value. The fitness 
value of each chromosome is assessed and a survival of 
the fittest strategy is applied. There are three main 
genetic operators of the GA, these are known as Selection 
or reproduction, crossover and mutation. 

The steps involved in implementing the GA are as 
follows: 

1. Generate an initial, random population of 
individuals for a fixed size. 

2 . Applying it to the process 

3 . Evaluate the process using fitness function chosen 
i.e. IAE plus weighted control efforts (squared) 

4 . Select the fittest members of the population. 

5 . Reproduce using a probabilistic method. 

6. Implement crossover operation on the reproduced 
chromosomes 

7 . Execute mutation operation with low probability. 

8. Repeat step 2 until a predefined convergence 
criterion is met. The convergence criterion of a 
genetic algorithm is a user- specified condition e.g. 
the maximum number of generations or when the 
string fitness value exceeds a certain threshold. 

An illustrative flowchart of the GA implementation is 
shown in Figure 1 

3. Genetic Operator 

To maintain the genetic diversity GA uses the operators. 
Genetic diversity is a necessarily procedure for process of 
evolution. Commonly three main Genetic operators of 
selection or reproduction, crossover and mutation are 
used in implementation of GA. 

A. Selection or Reproduction 

Selection is usually the first operator applied on 
population. The Selection operator is also known as 




reproduction operator. Selection means extract a subset 
of genes from an existing population. Every gene has a 
meaning so one may extract a kind of quality 
measurement called as fitness value and by following this 
quality selection can be done. Fitness function quantifies 
the optimality of a solution so that a particular solution 
may be ranked against all the other solutions. The 
function shows the degree of the closeness of the given 
solution to the desired solution. Roulette wheel selection, 
Rank selection, Boltzmann selection, Steady state 
selection and tournament selection are examples of 
selection operators [18]. 

B. Crossover 

After the selection, crossover operator takes place and it 
combines two chromosomes (parents) to produce anew 
chromosome (off spring). The centre idea of crossover is 
that the new generated chromosome may be better than 




Figure 1 Flowchart of Genetic Algorithm 



both the parents if it takes the better characteristics from 
each of them. It depends on user specified crossover 
probability. It indicates how often crossover is performed. 
A probability of 0% means that the ‘offspring’ will be 
exact replicas of their ‘parents’ and a probability of 100% 
means that each generation will be composed of entirely 
new offspring. One point crossover, two points, uniform, 
arithmetic and heuristic are few types of crossover. 

C. Mutation 

Mutation is the last operation in GA. It is used to 
maintain genetic diversity from one generation of a 
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population of chromosome to the next generation. It 
depends upon the user defined mutation probability. 
Mutation probability values of around 0.1% or 0.01% are 
common, these values represent the probability that a 
certain string will be selected for mutation i.e. for a 
probability of 0.1%; one string in one thousand will be 
selected for mutation. Mutation alters one or more gene 
values in a chromosome from its initial state. This can 
result a fairly new gene value being added to the 
population and with new gene value the GA is able to get 
better solution in contrast to previous solution. Mutation 
helps to prevent the population from stagnating at any 
local minima. The mutation operators are of many types 
such as flit bit, Boundary, Non uniform, Uniform and 
Gaussian. 

4. Process Evaluation Based on Fitness 
Function 

A. Fitness function / Performance index 

The most crucial action to choose the fitness function to 
evaluate the fitness of each chromosome. Performance 
indices such as mean square error (MSE), integral time 
multiplied by absolute error (ITAE), integral absolute 
error (IAE) and integral square error are commonly used 
as the fitness function in applying GA [7]. 

B. Fitness Value 

The PID controller is designed to minimize the fitness 
function which is function of system error. The smaller 
value of fitness function indicates the higher fitness value 
of the corresponding chromosome and vice versa, and it 
is represented by 

fitness value = (1) 

fitness function 



u(k)= k p (k)Xj + k i (k)x 2 + k d (k)x 3 (4) 

n 

where x 3 (k) — e(k) , x 2 = ^ e(k)T , 

k=l 

e( k) - e( k - 1 ) 

x 3 = — , and k p , k h k d are proportional, integral, 

and derivative gains respectively and T is sampling time. 
The convergence criterion of a GA is a user- specified 
condition and termination of tuning of PID parameters 
can be done by reaching the maximum number of 
generations or when the fitness value exceeds a certain 
threshold. To get the satisfactory dynamic characteristics 
IAE has been considered as the fitness function in the 
present work. To limit the controlled input from being 
too large, it also accumulates the quadratic term of 
controlled input to the fitness function. The expression 
for modified fitness function can be represented as 
follows 

J = Y J (w 1 \ e (k)\ + w 2 u 2 (k)) (5) 

k=l 




5. Problem Formulation 

Figure 2 depicts the schematic diagram of the proposed 
genetic tuning of PID parameters. The scheme composes 
of usual feedback control system which has the full 
model of coupled nonlinear multivariable plant and the 
PID controller. The parameters of controller have been 
recursively tuned online by GA. To begin with, we 
assume that system to be controlled has multiple inputs 

u(k) E y{ mX 1 and multiple outputs y(k ) E tyi p x 1 and 
the control objective is to follow a prescribe set points 
r(k) E x 1 . Let the tracking error be defined by 



e(k) =r(k)-y(k) (2) 

and the controlled input of the system is 

u(k) = K(k)X(k) (3) 



where K(k)E is the set of controller’s 

parameters \k p (k) k t (k) k d (kj] 



and 



of and 



X(k)eW xp is the set of errors which is represented 
by \x,(k) x 2 (k) x 3 (k)f. 



Figure 2 Schematic diagram of PID tuning based on GA 



6. The Turbo - Generator Plant 

The schematic diagram of laboratory scale turbo- 
generator plant is shown in Figure 3, It was developed as 
an emergency power supply for biological experiments 
[6]. A compressed centre air supply runs a turbine which 
drives a Y - connected three phase synchronous generator. 
The field of the generator is separately excited by an 
external source. A balanced three phase resistance serves 
as a load which can be switched between 1.5 Q and 2.5 Q. 
A tacho generator is used for speed measurements. The 
percentage of valve opening V (%) and field current i E 
(mA) are input signals Ui(t) and u 2 (t) respectively. The 
turbine speed n (min' 1 ) and the rectified mean value of 
the generated voltage U A (V) are considered as the output 
signals yi(t) and y 2 (t) respectively. The numbers of 
machines poles are constant so the frequency f (Hz) of 
the generated unrectified voltage is proportional to speed 
n. 
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Figure 3. Schematic diagram of laboratory scale turbo-generator plant 



different techniques is also shown. It can be seen that 
with the help of proposed scheme system response 
becomes quite faster, rise time becomes shorter and 
smaller controlled input is required to meet the desire 
response (Figure 5). The optimization process of the 
fitness function is shown in Figure 6 which shows the 
optimize value of fitness function at each iteration. 
Finally optimize PID parameters 

" 0.1985 0.1458 0 . 0119 " 



are K = 



0.1166 0.0454 0.0207 



All the computation has been done on Intel Core 2, 

| 3 gp@ 1.86GHz machine and a fair comparison of 

computational burden is shown in Table 2 which 
indicates that due to complexity, proposed algorithm 
consumes more time to tune the PID parameters. 



7. Simulation Results 

Turbo- Generator plant which is nonlinear coupled 
MIMO system is used to illustrate the effectiveness of the 
proposed scheme. The Input-Output relation is as follows: 

=0.4y,(k-l)+ Ul ^ 1} +0.2u l \k-l)+0.5u 2 (k-l) (6) 

1+Uj (k-1) 

y 2 (k) = 0. 2y 2 (k- 1 )+ +0.4 m 2 3 (£-1)+0.2w 1 (£-1) (7) 

X 1 14 2 yK” ly 

The GA is used to tune for the optimal PID parameters 
that minimize the fitness function (Eq. 5). Therefore the 
parameter tuning task of the PID controller using GA can 
be considered by selecting the three parameters k p , ki and 
k d such that the response of plant will be as closest to 
desire response. The GA parameters chosen for 
simulation are shown in Table 1. 



Table 1 Parameters of GA 



GA property 


Value/ Method 


Population Size 


30 


Maximum Number of 
generations 


100 


Fitness Function 


Quadratic form as in Eq (5) 


Selection Method 


Roulette wheel selection 


Crossover Probability 


0.90 


Crossover Method 


Arithmetic Crossover 


Mutation Probability 


0.01 



Initially many individual solutions are randomly 
generated to form an initial population, allowing the 
entire range of feasible solutions that results in stable 
response (0<k p i<0.7, 0<kii<0.9, 0<k dl <0.3, 0<k p2 <l, 
0<k i2 <l, 0<k d2 <l). Initial solution can be generated by 
following expression 



x = a+ (b-a)*rand 



( 8 ) 



where a, b are upper and lower limits of the parameter 
value k p , ki and k d and rand is a random number between 
0 to 1 . System response to sum of delayed step function 
is shown in Figure 4 and a fair comparison among three 



Table2 Computational burden 



Technique 


Computational time( sec) 


Fixed PID 


0.0290 


DRNN based PID Tuning 


0.0500 


GA based PID Tuning 


3.6020 




Figure 4 System response to sum of delayed step function 




contrafod input ul by GA based tuned PID 
ronliotod input ul by DRNN bated tuned PID 

condoned input ul by fixed P© parameters 

runlinOrd input u? by GA bated tuned PID 

cmrolled input u2 by 0RNN based tuned PC 
controlled input u2 by Fixed P© parameters 



Figure 5 Controlled inputs 
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8. Conclusions 

This paper demonstrates that how GA can be used to 
optimize the parameters of PID controller for turbo- 
generator plant. The features of the proposed scheme are 
summarized as follows: 

• Optimized parameter can be obtained by optimizing 
the PID parameters controller with GA. 

• It is important to select permissible upper and lower 
limits of parameters otherwise an optimal solution 



is not possible in realistic time frame i.e. online / 
real time. 

• The parameter optimization is global and has there 
is no connection with initial population. 

• With lower rise time and controlled input the 
proposed scheme is capable to better rack the 
delayed step function compare to other techniques. 

• Controlling input effort is reduced and this offers a 
major advantage. 
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Abstract 

In this paper an observer-based stochastic controller is 
developed for linear stochastic discrete-time dynamical 
systems with uncertain parameters. A state estimator, 
based on Kalman filter, is firstly developed which 
incorporates in its structure the statistical information of 
the uncertain parameters. Although the developed 
estimator uses the nominal values of the system 
parameters; it showed to be stable, robust and gives 
satisfactory results. Such an estimator is then used to 
design an observer-based controller for discrete-time 
linear stochastic dynamical systems. The effectiveness of 
the theoretical derivation of the developed estimator and 
the observer-based controller are illustrated via a 
numerical example and a practical case study. 

Keywords: Discrete-time systems, linear systems, 

nonlinear systems, state estimation, stochastic systems, 
control design. 

1. Introduction 

Since most, if not all, real control systems are stochastic 
by their nature, it has been essential to develop and 
enhance stochastic control theory which mainly deals 
with this type of problems. With the progress of 
computational tools, and the development of state space 
approach, it became a systematic procedure to represent 
system models by state equations rather than frequency 
domain models. In practice, the states of the system are 
not accessible to direct measurements. Therefore, it was 
mandatory to develop a theory for state reconstruct of the 
system in order to generate the desired control strategies. 
For deterministic linear finite-dimensional time invariant 
systems, an observer was first introduced by Luenberger 
[1] leading to the asymptotic estimation of the state. In 
case of linear dynamical system with random 
disturbances, where the stochastic phenomena appears, 
the state estimation problem is solved using Kalman filter 
(KF) [2]. In standard Kalman filter, all the system 
characteristics (i.e. system model, initial conditions, and 
noise characteristics) have to be specified a priori. 
However, if there is an uncertainty in any of these 
characteristics, the filter may not be robust enough. For 



state estimation of linear discrete systems with uncertain 
parameters, interested readers can refer to [3] and 
references therein. On the other hand, for nonlinear 
dynamical systems, the stability, and hence the 
robustness of the estimator become worse. Several 
algorithms have been proposed to estimate the states of 
nonlinear systems [4, 5, 6, 7, and 8] and references 
therein. In general, the extended Kalman filter (EKF) is 
the most extensively used algorithm for state estimation 
of nonlinear systems. However, it is not the optimal 
estimator, and it may diverge if the system acquires 
strong nonlinearities, or the initial estimate of the state is 
wrong, or the process is modeled incorrectly, owing to 
linearization [9, 10, 11]. 

In this paper, a new approach is proposed to estimate the 
states of linear stochastic discrete-time dynamical system 
with uncertain parameters. The system model and the 
measurements are assumed to be corrupted by 
uncorrelated zero mean white Gaussian noise sequences. 
The parameters of the system are assumed to be uncertain, 
which obviously leads to an estimation problem of 
nonlinear (bilinear) stochastic discrete-time dynamical 
system. The new filter handles this problem and its 
performance is compared with that of the KF, when the 
parameters are assumed to be certain, and with the EKF, 
when the parameters are decided to be estimated. The 
proposed estimation technique is used to design an 
observer-based controller for linear stochastic discrete- 
time system with uncertain parameters. For linear 
stochastic systems with certain parameters, the separation 
principle holds. Therefore, controller design and state 
estimator can be performed separately. This leads to an 
optimal closed loop stochastic controller [12]. However, 
in case of stochastic systems with uncertain parameters, 
the separation principle does not hold, since additional 
terms will appear in the state space representation of the 
augmented system due to the uncertain components of 
the system parameters. Although this is the case, we can 
still design a closed loop suboptimal observer-based 
controller provided that it leads to an asymptotically 
stable system. To show the effectiveness of the 
developed estimation approach, illustrative examples are 
solved using the proposed estimator, the KF and the EKF 
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with parameter estimation. MATLAB is used as a 
programming tool and Monte Carlo simulation is applied 
to compare the behavior of the different filters. An 
observer-based controller, using the developed estimator, 
is designed to regulate a DC motor. The obtained results 
show the applicability and efficiency of the proposed 
procedure in handling this problem. 

The rest of the paper is organized as follows. In section 2, 
the control problem for stochastic linear discrete-time 
dynamical systems with uncertain parameters is 
formulated. Section 3 is devoted to the presentation of the 
mathematical structure of the developed estimator. The 
dynamics of the observer-based controller is presented in 
section 4. The algorithm used to implement the closed 
loop stochastic controller is given in section 5. In section 
6, illustrative examples are solved to show the 
effectiveness of the proposed design approach. The paper 
is concluded in section 7. 



2. Problem Formulation 

Let us define the following stochastic control problem: 



min J = E ■ 

Subject to: 
x 



1 kf ~\ 
k=k„ L 



WXu -x 



II M 



+ H, 



* + 1 ^ A k*k + B k“k +d +w k 



d-a) 

d-b) 

(1-c) 

is the 



yjc + 1 - Cjc+l x k+l + v k+l 

where: x k g R n is the state vector, y k eR n 
output vector, u k eR r is the control vector, d g R n is 
a constant vector, x d eR n is the desired state vector, 



Au — 



M 1 



j^nxn 

a- 



is the system matrix with uncertain 
parameters, ;/,y g {l,2,....,n} assumed to be 

uncorrelated white Gaussian random variables with the 
following statistics a t j k -N(da,<j^ ), B h e R nxr is the 



u 



control transition matrix, w k g /? " is a zero mean 

white Gaussian input noise vector with covariance matrix 
Qk = E{w k wl } e R nxn , C k+l = [c ijk+i ] e R mxn is the 

output measurement matrix with uncertain 
parameters, c ij k+i J e {1,2 m},y g {1,2 ,...., w} assumed 

to be uncorrelated white Gaussian random variables with 
the following statistics 
a zero mean white Gaussian output noise vector with 
covariance 
M g R nxn 

the states, N g R rxr is positive definite weighting matrix 
for the control, x Q g R n is the initial conditions of the 
states assumed zero mean random Gaussian vector with 



■N(c ij ,C7l),v k+l sR m is 



matrix R k+l = E{v k+l v k+l } e R mxm 
is positive semi-definite weighting matrix for 



R n 



and 



covariance matrix P 0 =P 0 \ 0 =E{ x 0 x 0 

finally k g {0,1,2,...} is the discrete time instant. 

2 

Let b k g R n ,c k+l g R nm be the parameter vectors 
defined as: 



\ ~ ’ a 12 k ’-’ a 21 k > a 22 k 

C k + 1 = [ C H t+ i ,C 12 i+1 » — > C 21 i+1 > C 22 t+1 >-f 

In the above model it is assumed that 
x o>bk’ c k+ 1 > w k’ v k+ i a re all independent. Moreover, it is 
assumed that the model has the following properties: 
E{x k w ] } = 0, V* < j ; E{x k v ] } = 0, Vft, j ; 

E{y k v] } = 0, V* < j ; E{x k b ] } = 0, V* < j ; 

(3) 

E{x k c] } = 0, ; E{y k c ] } = 0, Vft < j; 

E{b k c]) = (),VkJ 

Since the parameters of the system are uncertain, then 
we have a bilinear estimation problem rather than linear, 
which intern leads to a non-Gaussian estimation problem. 
Our objective is to get a good estimator for the state 
vector, and hence use the estimated states to generate the 
desired control signal. 



3. The Developed Estimator 

Let us for the moment consider the following problem: 

**+i = A kX k +W k (I'-a) 

y k+ \ =C k+l x k+l +v k+1 (l'-b) 

The system model (I') has the same definitions and 
properties of the system model (1) in section 2. 

State estimation problems with uncertain parameters is 
usually solved by treating the parameters as static states, 
i.e. P k+l = p k with the prior at k - 0 being N(0,Z) , 
where /? is the vector of the mean or nominal values of 
the parameters and I is its covariance matrix. This 
approach increases the dimensionality of the problem, 
and hence the computational burden. To avoid this 
difficulty, the nominal values of the parameters will be 
used in the estimator. However, the impact of this 
approximation on the dynamics of the estimator will be 
investigated. 

Assume that the filtered estimate x m and its associated 

covariance matrix P k \ k are given at the k th sampling 
instant, and used to approximate the conditional 
probability density function p(x k I Y k ) by Gaussian 
distribution with conditional mean and covariance matrix 
as given above. 

The filtered estimate of system (I') and its associated 
covariance matrix are given by: 

a) The predicted estimate of the state vector x k+Vik , its 
associated covariance matrix P k+ \\ k , and the predicted 
estimate of the output vector y k+ \\ k , are given hy: 

4+n k = A h\k ( 4 ) 

4+11* = AP k\k AT + 4 1* 4 k ®k\k + -4 + Qk (5) 
.y*+ii* = Cx k+m (6) 

b) The filtered estimate of the state vector x k +\\k +\ , its 
associated covariance matrix P k+ i\ k+ i , and the filtered 
estimate of the output vector y k+ \\ k+ \ , are such that: 
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x k+m+i = x k+m + K k+iyjc 



+11 k 



(7) 



K k+ 1 ~ P k+l\k C T [CP k+llk C T + ^k+m^k+i^k+wk + 1 + ^k+ 1 ] 1 

(8) 

(9) 



( 10 ) 



^A:+1IA:+1 = ~~ ^fc+l^^fc+llA: 

JVA:+1IA:+1 = C^fc+llfc+l 
Proof: 

a) Derivation of the Prediction Estimator and the 
Corresponding Covariance Matrix: 

Equation ( 1 -a) can be rewritten as: 

x k + 1 = + Ak )(*k\k +x m) + w k 

x k + 1 = Ax k\k + Ax k\k + A k *k\k + + w k 

where A is the matrix contains the nominal (mean) 
values of the parameters, while A k contains the 
perturbations around these nominal values. 

/«*/ rwi 2 

Let b k = \a n ,d 12 >&22 k *— T e be a vector 

contains the deviation of the parameters of the matrices 
A k around their nominal values A , and define: 

O O 



x k\k 



®k\k ~ 



®k\k ~ 




gR u 



( 12 ) 



gR n 



where O = [ 0 0 0 is a vector of zero 

elements. 

Then, (11) can be rearranged in the following form: 

**+i = A**i* + Ax m + S m b k + 0 k]k b k +w k (13) 

The predicted value of the state vector is given by: 



= E{x k+1 \Y k 



*k+i\k ~ ^i^k+i 

Using the model properties (3), the predicted estimator is 

such that: 

*k+\\k = Ax k \ k (15) 

From (13), (15), the prediction error is given by: 
*k+\\k = Ax k \ k + 0 k \ k b k + 0 k \ k b k + w k (16) 

The covariance matrix of the estimation error is such that: 

P k+i\k = ^{**+ 11 ** 1 + 11 *} ( 17 ) 

Using the model properties (3), p k +\\k is given by: 

(18) 



P k+l\k ~ AP k\k AT + Kk E A +S k+Qk 



r k =diag[al t a Un <x ( 

S k =diag[e Ut 

n 

e " = 2 P k\k„ 



12 * 

e 22 , 



where: 

(19-a) 

(19-b) 

(19-c) 



Equation ( 1 -b) can also be written in the form: 

yk + 1 =(C + Ck+0( x k+l\k + x k+l\k ) + v k + 1 

yk+l = C x k+l\k +C x k+l\k + Ck+\ x k+l\k +C k+l x k+l\k +v k + 1 

( 2 °) 

where C is the matrix contains the nominal (mean) 
values of the parameters, while C k+l contains the 
perturbations around these nominal values. 

Let Cu, i =[cii ,c 19 ,...,c 91 ,..] r g R nm be the vector 

lt-rl L1 k + 1 AZ '*+ 1 Zji k+\ 

contains the deviation of the parameters of the matrix 
C k+1 around their nominal values C . Then, by 
rearranging (20), one gets: 

yk+l = Cx k+l\k + Cx k+l\k + ^k+Wk^k+l + ^k+Wk^k+l +V A:+ 1 

( 21 ) 



where: 



5Pj 



k+\\k ~ 



v k+\\k 



o 



o 

x k+\\k 



o. 

o. 



o o 



" x k+l\k 



gR 



<mx(nm ) 



r, 



k+l\k 



x k+l\k 



o 



o 



ZT 

x k+l\k 



o. 

o. 



( 22 ) 



R 



nx(nm) 



O O x k+i\k 

From (21), the predicted output is given by: 

yk+l\k = C x k+l\k 

Accordingly, y k+ \\ k is such that: 

yk+l\k = C x k+l\k +l ^k+l\k^k+l + ^k+l\k^k+\ +v k + 1 

Hence: 



Jjt+iifcJyt+iut 



- Eih+Ukh+Uk 1 



- EP k+\\k E + 'f'k+Wk^ k+^k+Wk + E k+\ + R k + 1 



(23) 



(24) 



(25) 



where: 



(14) U k+ 1 =<Bag[al a ] 



(26) 



^12fc + i 

L k+ \ =diag[l lh Z 22ti Wl 

** + , =±Pk 

7=1 

b) Derivation of the Filtered Estimator and the 
Corresponding Covariance Matrix: 

The filtered estimate of x k+l is given by: 

x k+i\k+i = E{ x k+i ,y k+ il 



= x k+l\k + r> 



>-i 






yk+l\k 






= E {x k+ \y k+ \\k) 



(27) 



(28) 



= E Vx k+ \\k + Xk + \\k)h + \\k} 

Using (24) in (28), and due to the independency 
between x k+M ,c k+l ,v k+1 , one gets: 
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Let 



K k+ 1 



= p - pr 1 _ 

•* 7 + 1 37 + 11 * 37 + 11 * 37 + 11 * 



(30) 



Using (25), (29) in (30) one can get K k+1 as given by 
(8). Hence: 

X A;+1IA;+1 = x k+l\k + K k+l h+llk (31) 

The filter estimation error is such that: 

x k+ i\k+i = x k+m ~ K k+ ih+m (32) 

Using (32), the covariance matrix of the filter estimation 
error is such that: 

B k+ m+i ~ Kk+iC)Pk+i\k (33) 



4. Control Design 

In the system model (1), since parameters of the system 
as well as the states are random, the system is bilinear 
rather than linear. Therefore, the separation principle 
does not hold. However, the problem will be treated as 
deterministic and the resulted control strategy will be 
applied to the stochastic problem. It is obvious in this 
case that the resulting closed loop control system is 
suboptimal rather than optimal. 

Let us consider the following servo-mechanism discrete- 
time linear deterministic system: 




Subject to: 

x k + i =^ x k + B k u k + d (35) 

where x k eR n is the state vector assumed to be 

deterministic, A e R nxn is the matrix contains the 
nominal values of the parameters, while other variables 
are as defined in section 2. 



By writing the Hamiltonian of this problem, and from the 
necessary conditions of optimality, the expression for the 
control is such that [13]: 



u k = 


-N~ 


l B T k A 




{p rk 


-M)x k+ Z k+ Mx d ] 


(36) 


where P r 

r k 


= M 




p, 

r k 


3 “‘A 

+1 


with P r = [0] 

k f 


(37) 


with 












3= / +B k N~ l B 


T P 

k r k+l 


(38) 
















4 = 


-A 7- 


P r 3 

r k+l 




N~ 


‘^4+1 ~Mx d +A r P rku 


3 ~ x d 






<vk+ 1 






with 


o 

II 


(39) 


Since 


the 


state 


vector 


x k is 


not available for 


direct 



measurement, the estimate of this vector x klk will be 

used in conjunction with the state feedback controller 
(36). Therefore, the dynamics of the observer-based 

closed loop stochastic controller at the k th iteration can 
be summarized as follows: 

a) Calculate, off-line, P rk , for k = k 0 ,....,kj_ x using 
(37)-(39). 

b) Calculate the control signal using (36). 

c) Calculate the predicted value of the state vector as 
given by: 

=*7i* + B k u k (40) 



d) The covariance matrix P k+i \ k , and the predicted 

estimate of the output/ k +m are as given by (5), (6). 

e) The filtered estimate of the state vector x k+m +1 , its 

associated covariance matrix P k+m+x , and the 
filtered estimate of the output vector y k +m +1 are as 
given by (7)-(10). 

The flowchart of the observer-based controller is 
given in the following section. 



5. The Flowchart of the Algorithm 




No 



END 





6. Simulation Results 

6.1 Illustrative Example : 

To investigate the effectiveness of the developed 
estimator, we employ firstly the following numerical 
values of a three dimension (3-D) discrete-time linear 
system model: 



A = 



1 

-0.009889 



-1.98 



0.009983 

0.9949 

-0.9874 



0 

0.0091 

0.8305 



,C = [ 1 0 0], 



X 0 =[l 1 lf,x olo =[0 0 of. 



P o\o - 1 where I e R 3x 3 is the unity matrix. 

The states of the system are estimated using the 
developed filter, KF, and EKF. Different values of 
Q , R and different levels of uncertainties in the system 

parameters are used in our simulation. For comparison 
purposes, Monte Carlo simulation is applied to explore 
the qualitative behavior of the root mean square errors 
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(RMSE) of the estimated states resulting from the 
application of the three filters. This is defined by: 



RMSE {xj (k ),Xj (k \k)) = 



(x j (k ] ~ x j (k lk )}2 



MI 



for \/k = f and j = l,....,n (41) 

where MI is the number of Monte Carlo iterations. 
Moreover, the root mean square error index (RMSE 
index) for each state is calculated using the following 
formula: 



RMSE index (xj ) = 



\ ( k f )(MI ) 



(42) 



Table (1) shows the RMSE indices calculated over 100 
Monte Carlo iterations for each case study while using 
the three filters. 



Table 1: RMSE index (MI = 100) 



Case 


X) 


Q 


R 


Pn% 


RMSE index 












MKF 


KF 


EKF 


a 




0.1 


0.1 


10 


0.2597 


0.269 




b 


1 


0.5 


1.0 


20 


0.7947 


1.821 


Div. 


c 


0.5 


1.0 


30 


0.8235 


8.233 


d 




1.0 


1.0 


30 


0.8717 


26.444 




a 




0.1 


0.1 


10 


1.8834 


1.887 




b 


2 


0.5 


1.0 


20 


19.687 


25.099 


Div. 


c 


0.5 


1.0 


30 


61.726 


66.207 


d 




1.0 


1.0 


30 


576.66 


686.21 




a 




0.1 


0.1 


10 


10.054 


10.09 




b 


3 


0.5 


1.0 


20 


95.214 


125.57 


Div. 


c 


0.5 


1.0 


30 


299.092 


345.97 


d 




1.0 


1.0 


30 


2.23e+3 


2.87e+3 





From Table 1, the proposed filter has the best 
performance in all cases. It has the minimum RMSE 
index calculated over 100 Monte Carlo iterations, 
whereas the (EKF) fails to estimate the states since it 
diverges. The RMSE plots of the states using the three 
filters for Case-a are shown in Figures (1-3), and for 
Case-b are shown in Figures (4-6). 



Modified Kalman Filter 



Kalman Filter 




0 1 1 I I I I I I I I I 

0 1 00 200 300 400 500 600 700 800 900 1 000 

time (sec) 




Figure 2: RMSE ofx2 (Case-a). 



Modified Kalman Filter 




Figure 3: RMSE of x3 (Case-a). 



Modified Kalman Filter 



t r ^ ~ n ; r^r ^ i ;«i ~ i i ”i 




0 100 200 300 400 500 600 700 800 900 1000 



Kalman Filter 



T 1 I 1 I “ T ~ 1 I ~ 1 




Extended Kalman Filter 



o 1.2 -4 - 4* t — -| — — — +- — — — I— — 1 - - 1 1 

s o 

i i i i i i i i i i 

0 100 200 300 400 500 600 700 800 900 1000 

time (sec) 



Modified Kalman Filter 




Figure 1: RMSE of xl (Case-a). 



Figure 4: RMSE of xl (Case-b). 



Modified Kalman Filter 

200 




0 I I I 1 I I L-l I I I I 

O 100 200 300 400 500 600 700 800 900 1000 

time (sec) 



Figure 5: RMSE of x2 (Case-b). 
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Modified Kalman Filter 



600 




Kalman Filter 




0 1 i I I I LJ I I I I 

0 100 200 300 400 500 600 700 800 900 1000 

time (sec) 



Figure 6: RMSE of x 3 (case-b) 



6.2 Case Study: 

The transfer function of a DC motor is given by [14]: 



K r 



m _ 

E(s) s [(sL a + R)(sJ + B) + K T K b ] 



(43) 



where, 

6 is the motor angle displacement. 

6 is the motor rotating speed, where 6 = — . 

dt 

i a is the armature winding current. 

E is the applied voltage. 

K t is the motor torque constant. 

K b is the back emf constant. 

L a is the armature winding inductance. 



R is the armature winding resistance. 

J is the moment of inertia of rotor and load. 

B is the damping coefficient. 

The state vector x is defined by x T -[6 6 i a f . 



For the above model, it is assumed that the rotating 
speed of the motor 0 is measured at the output. 

6. 2. a) State estimation under no control: 

The continuous state equation of the DC motor is: 



d_ 

dt 



9 ' 




"0 


1 


0 


' 9 ' 




" 0 




6 


= 


0 


-B/J 


k t /j 


9 


+ 


0 


E (44) 


ja _ 




0 


~K b /L a 


~E/L a _ 


Ja _ 




4 





y = [0 l 



0] 



6 

0 



(45) 



where £=150 is the reference input. 

The above model is discretized using Euler’ s method 
with sampling period Ts= 0.005. Therefore, the discrete 
model is such as: 



4+i 




"1 


0.005 


0 


4 




0 


4+i 


= 


0 


0.9994 


0.0749 




+ 


0 


J'ak+l _ 




0 


-0.9632 


-0.9238 


jak _ 




0.5023 













150 (46) 



yk+ 1 - to i 



0] 



4+i 

4+i 



W i 



(47) 



where, the system model is corrupted by an input noise 
w (k ) ~ N (0, 0. 1) and measurement disturbance 



v (Jfc) ~N (0,0.01). 



The system parameters K T ,K b , L a and R , are 
assumed to be uncertain with a standard deviation equals 
to certain percentage P n of their nominal values. 

x 0 = [10 80 0.5f+ olo =[5 79 1 ] T , P olo = I e R 3x3 

The RMSE indices of the states calculated over 100 
Monte Carlo iterations using the three filters are shown in 
Table 2. RMSE plots for Case-b are shown in Figures (7- 
9). 



Table 2: RMSE index, Q=0.1, R=0.01 (MI = 100). 



Case 


State 


Pn% 


RMSE index 








MKF 


KF 


EKF 


a 


9 


5 


8.6057 


8.6057 




b 


10 


9.3070 


9.3070 




a 


9 


5 


0.1034 


0.1231 


Div. 


b 


10 


0.1044 


0.2213 


a 


ia 


5 


6.6085 


12.4121 




b 




10 


15.5597 


34.4408 





Modified Kalman Filter 





time (sec) 



Figure 7: RMSE of 6 (case-b) 
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Modified Kalman Filter 



100 200 300 400 500 



700 800 900 1000 



Kalman Filter 




0 I ' 1 1 1 1 1 1 1 

0 100 200 300 400 500 600 700 800 900 1000 



time (sec) 



Figure 8: RMSE of & (Case-b). 



Modified Kalman Filter 




0 | 1 1 1 1 1 1 1 1 1 

O 100 200 300 400 500 600 700 800 900 1000 

time (sec) 



Figure 9: RMSE of i a (Case-b). 



From the simulation results, the proposed estimation 
procedure showed to be the best, and leads to the 
minimum value of the RMSE index. Therefore, it will be 
used in the implementation of the observer-based 
controller as shown in Figure (10). 




Figure 10: Block diagram of the closed loop system. 



6.2.b) Observer-based Controller. 

For the problem in hand, it is desired to regulate the 

system to the rotating speed 6 d =76.8085 rpm . The 
initial values for the states, the estimated states, and 
covariance matrix P are as follows: 

x 0 =[5 40 lf,x olo =[() 0 0 f,P oio =I eR 3x3 

The state weighting matrix is M - diag [ 0 1 1] , the 
control weighting matrix is N =[1], the covariance matrix 
of the input noise is 2 = 0.17, the covariance matrix of 

the output noise is/? = 0.01, and the percentage change 
in the parameters is P n = 10% . The actual and estimated 

closed loop trajectories as well as the control signal are 
shown in Figures (11-14). 







Figure 1 1 : Actual and Estimated motor angle 
displacement 0. 




Figure 12: Actual and Estimated motor 
rotating speed ^ . 



!■ " ' 




am ■ xheO 



Figure 13: Actual and Estimated armature 
winding current^. 



I.‘ : :r*: <: mar F lier 




4CC SB 600 
:ine k: 



Figure 14: Feedback Control u. 
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7. Conclusion 

In this paper, an observer-based controller is developed 
for linear stochastic discrete-time dynamical systems 
with uncertain parameters. A state estimator, based on 
Kalman filter, is firstly developed which incorporates the 
statistical information of the uncertain parameters in its 
structure. Although this problem is treated in the 
literature as a nonlinear estimation problem; with the 
developed approach the system is still treated as linear 
and without any increase in its dimensionality. As a result, 
we have less computational burden and better numerical 
properties. Based on the developed filter, an observer- 
based feedback controller is proposed. Simulation of 
different illustrative examples showed that the state 
estimation using the proposed filter is better than KF, 
while the EKF diverged and failed to estimate the system 
states. Moreover, although the separation principle does 
not hold, simulation results showed that the regulated 
system is asymptotically stable, and the states of the DC 
motor reached the desired ones. However, it is still 
necessary to analyze the stability of the controlled system 
more rigorously in a theoretical frame work. This 
problem will be treated in the future. 
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Abstract 

In this paper a novel nine-level inverter system for Direct 
Torque Control (DTC) induction motor drive is presented. 
DTC strategy of a nine-level Multi Point Clamped (MPC) 
Voltage Source Inverter (VSI) fed induction motor is 
implemented. Nine-level MPC inverter is an advanced 
technology and are commonly applied to high power 
electrical drives instead of the two-level VSI.Nine-level 
inverter has 729 space vector switching states and there 
are 106 effective vectors are possible. The proposed nine- 
level inverter scheme is capable for enough degrees of 
freedom to control both electromagnetic torque and stator 
flux with very low ripple. From the simulation results 
shows that feeding electrical drive with nine-level 
inverter can greatly improves the drive performance as 
compared to the seven-level inverter. The performance of 
this control method has been demonstrated by 
simulations performed using a versatile simulation 
package, Matlab/Simulink. 

Keywords: DTC, Nine-level inverter Multi point 

Clamping, Switching table, Induction motor, 
Matlab/Simulink. 

1. Introduction 

Multilevel converters have recently increased interest in 
the research and industry communities [3]-[7]. In these 
kind of converters, the output voltage can take several 
discrete levels of equal magnitude. The multilevel 
converter first proposed in [1] was aimed at reducing the 
harmonic content of generated voltage or current 
waveforms. The harmonic content of such a waveform is 
greatly reduced, if compared with a two-level waveform. 
Nowadays, there are three prominent multilevel 
topologies. Multilevel structures derived from those 
proposed in [2] and [3] are known as diode-clamped (or 
neutral-point-clamped) and flying-capacitor multilevel 
converters, respectively. Both of them are based on a 
serial connection of low- voltage power semiconductors. 
The use of voltage clamps to equally share the input 
voltage between the serial-connected devices and at the 
same time having a multilevel output voltage is a 
common feature of these two structures. Thus, by this 
means, the known problems posed by the serial 



connection of power semiconductors are avoided and a 
multilevel output voltage is also possible. 

Feeding electrical motors with multilevel inverters can 
greatly improve the motor performance. The harmonic 
content of the voltage applied to the motor’s terminals 
will be greatly reduced. Recently, the high dv/dt 
produced by modem power electronic based inverters has 
been identified as the source of several motor failures [8], 
[9]. The ability of multilevel inverters to deliver a 
multistepped voltage waveform will prevent motor 
failures by reducing the voltage change rate of the 
applied voltage. Any of the three prominent multilevel 
topologies can be used in motor drive applications. 

MPC is so called since in their architecture there are 
several points clamped to specific voltages using some 
components. Even diode-clamped converters belong to 
this family because the bus between two switches is 
clamped by a clamping diode. Furthermore, when the 
number of voltage level is odd, the converters are called 
Neutral Point Clamped (NPC) because the neutral point 
is clamped, has found wide application in medium- 
voltage high-power applications [10] [11]. The main 
features of the MPC converter include reduced dvldt and 
Total Harmonic Distortion (THD) in its AC output 
voltages in comparison to the conventional two level 
converters. As in any multilevel converter it can be used 
in the medium-voltage applications to reach a certain 
voltage level without series connection of power 
semiconductors. 

In principle, DTC method is based on instantaneous 
space vector theory. By optimal selection of the space 
voltage vectors in each sampling period, DTC achieves 
effective control of the electromagnetic torque and the 
stator flux on the basis of the errors between theirs 
references and estimated values. It is possible to directly 
control the inverter states through a switching table, in 
order to reduce the torque and flux errors within the 
desired bands limits [12] [13]. The present work is based 
on the study of the application of DTC to the nine-level 
MPC VSI, and the advantages that can be obtained from 
using this topology when compared to the seven-level 
inverter. 
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The paper is organized as follows: the concept of nine - 
level MPC VSI is discussed in section 2 and 
representation of proposed DTC scheme is discussed in 
section 3. The simulation results and discussion of the 
proposed method are exposed at the next section. Finally 
conclusions are given in the last section 

2. Nine-Level MPC VSI and Related Output 
Voltage Vectors 

The architecture of MPC is several points clamped to 
specific voltages using some components. Even diode- 
clamped converters belong to this family because the bus 
between two switches is clamped by a clamping diode. 
Furthermore, when the number of voltage level is odd, 
the converters are called Neutral Point Clamped (NPC) 
because the neutral point is clamped. As Figure. 1 shows, 
the 9-level leg is completely different: in MPCs the 
voltages are clamped using couples switch-diode instead 
of using a simple diode. Anyway, given a number of 
levels, the number of switches needed by MPC is the 
same needed by diode-clamped. The control of MPC leg 
is more complicated in the respect of other topologies. 
Even this kind of converter allows finding 
complementary couples of switches. The constrain so 
introduced is not physiologically necessary, but it is a 
simple way to simplify the control scheme and the 
switching table (Table I). 

TABLE I 

SWITCHING STATES OF A NINE-LEVEL 
INVERTER 



Tj 


t 2 


t 3 


t 4 


Ts 


t 6 


Ti 


t 2 


t 3 


t 4 


Ts 


T« 


V 


1 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


4Uo 


0 


1 


1 


1 


0 


1 


1 


0 


0 


0 


0 


0 


3Uo 


0 


0 


1 


0 


0 


1 


1 


1 


0 


0 


1 


0 


-2Uo 


0 


0 


1 


0 


0 


1 


0 


0 


1 


1 


1 


1 


-Uo 


0 


0 


0 


1 


1 


1 


1 


0 


0 


0 


0 


0 


0 


1 


1 


1 


0 


0 


0 


1 


1 


0 


0 


0 


1 


0 


0 


1 


1 


1 


1 


1 


0 


0 


1 


0 


0 


0 


Uo 


0 


1 


1 


0 


1 


1 


1 


0 


0 


0 


0 


0 


2U 0 


0 


0 


1 


1 


1 


1 


0 


0 


1 


1 


0 


0 


3U 0 


1 


1 


0 


0 


1 


1 


1 


0 


0 


0 


0 


0 


4U 0 



Furthermore, to avoid shortcut, T 5 and T 6 must be 
complementary controlled. The same must happens for 
T' 5 and T' 6 . Table I is a possible switching table for a 9- 
level MPC leg; there is an intra-phase redundancy only 
for the middle level. Anyway, it is better to have T 6 
turned on in order to limit the reverse voltage drop 
uponT 3 .Dependently on the switching table used, the 
maximum reverse voltage drop over the components 
changes and a preliminary analysis must be done to 
choose the suitable component. Moreover, the switches 
in the middle of the leg must carry twice the voltage of 
the others. 

The DC bus voltage of 9-level inverter is split into nine 
levels by using eight DC capacitors Ci C 2 ,C 3 ,C 4 ,C 5 ,C 6 ,C 7 
and C 8 as shown in figurel. The voltage across the 
switches is l/8 th of the DC bus voltage, switching losses 



are cut in l/8 th with reduced harmonics of output 
waveforms for the same switching frequency, and power 
rating is increases. 

Relatively to the 7-level inverter which is only capable to 
produce 343 space vector voltages, a 9-level inverter has 
9 3 =729 switching states are possible as shown in figure. 
If voltages of eight capacitors are equivalent, some 
switching vectors are overlapped and there are 106 
effective vectors. 

According to the magnitude of the voltage vectors, we 
divided them into thirteen different groups: 




Figure. 1 : Nine-level multi point Clamped Inverter leg 
TABLE II 

Magnitudes of the voltage vectors 



Group 


Magnitudes of the voltage vectors 


1 


[V 0 ] 


2 


[Vi, v 2 , v 3 , v 4 , V 5 , V 6 ] 


3 


[V44, V45, V 46 , v 47 , V 48 , V 49 ] 


4 


[Vi2 3 , V 124, V125, V126, V127, V128] 


5 


[Vi 3 5, Vi 37 , Vi 3 9,Vi4i,Vi4 3 , V 145]? 


6 


[V 180, V 181, V 182, V 183, V 184, V 185, V 186, 
V 188>V 190, V 192, V 194, V 195] 


7 


[V 250? V 251 ? V 252, V 253? V 254, V 255? V 256? V 258? 
V 260? V 262? V 264? V 266] 


8 


[V 3 o 3 , V 3 Q4, V 3 05,V 3 06, V 3 07,V 30 8] 


9 


[V 3 95, V 3 96? V 397 , V 3 98, V 399 , V400] 


10 


[V 488? V 489? V 490? V 491 ? V 492? V 493 , V 494? V 495, 
V 496? V 498? V 500? V 502? V 504? V 506] 


11 


[V 588? V 589? V 590? V 591? V 592? V 593? V 594? V 595, 
V 596? V 598, V 600 V 602, V 604, V 606] 


12 


[V 676, V 677, v 678, V 679, V 680, V 681 , V 682, V 683, 

V684,V685,V686] 


13 


[V686, V 68 8, V 6 90, V 692 , V694, V 6 96] 
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The space vector modulation diagram of 9-level inverter 
is divided into six sextants, and each sextant is divided 
into sixty-four triangles regions in order to show the 
vectors nearest to the reference voltage as shown in 
figure 2 and 3. 

The voltage vectors of 9-level inverter are divided into 
thirteen groups as shown in table II. The procedure of 
selection of voltage vectors is similar to that of 3 -level 
inverter is discussed in [8]. 





Figure. 3: Zoom of the area near to the reference 
Space vector diagram of 9-level inverter in 
Sextant-I 



3. Representation of Proposed DTC Scheme 

The basic requirements of DTC strategy is the estimation 
of the stator flux and electromagnetic torque. The flux 
evaluation can be carried out by different techniques 
depending on the whether the rotor angular speed or 
position which are measured or observed. For multilevel 
applications, the voltage model is usually employed. The 
stator flux can be evaluated by integrating from the stator 
voltage equation 

a -«,<■> a* 



Further to calculation of the components of the flux, the 
estimated torque is determined from the following 
equation (2) 



T e =P(hs®ds ~L® q s) 



( 2 ) 



As shown on figure, DTC scheme needs only the output 
voltages and currents of the inverter which feeds the 
induction motor, the instantaneous values of flux and 
torque in the machine are then calculated and the error 
can be gotten after compared with the referring torque 
and flux (T eref and y/ ) 




F hex -angle y 



Figure.4: Configuration of proposed DTC Scheme 
The flux control is made by two-level hysteresis 

comparator and C 0 defines the action wished on the 
behaviour of the field. Also, a high level performance 
torque is required and the torque control is controlled by 
a hysteresis comparator built with eight lower bounds and 

eight upper know bounds and C r defines the action 



wished on the behaviour of the torque [9] . The hysteresis 
blocks are designed as it is shown in figure 5. 

For the switching vector selection, it is necessary to know 
the angular sector in which the actual flux is located. The 
actual position of the stator flux can be determined by 
equation (3), from the orthogonal flux components: 



6 = tan 1 



O , 

V as J 



( 3 ) 



The flux position in the (X and /? Planes is quantified in 
six sectors S of 60° degrees staring with the first sector 
situated between -30° and 30°. 



o' 




t i 


l 


-1 1 


r J 


L El 






Figure. 5: Flux and Torque Hysteresis blocks 
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TABLE VI 



Proposed Switching Table 





S 


Co 


Cf 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


+1 


-8 


684 


694 


685 


696 


680 


686 


681 


688 


682 


690 


683 


692 


-7 


594 


604 


595 


606 


590 


596 


591 


598 


592 


600 


593 


602 


-6 


494 


504 


495 


506 


490 


496 


491 


498 


492 


500 


493 


502 


-5 


254 


264 


255 


266 


250 


256 


251 


258 


252 


260 


253 


262 


-4 


184 


194 


185 


196 


180 


186 


181 


188 


182 


190 


183 


192 


-3 


127 


143 


128 


145 


123 


135 


124 


137 


125 


139 


126 


141 


-2 


307 


399 


308 


400 


303 


395 


304 


396 


315 


397 


306 


398 


-1 


5 


48 


6 


49 


1 


44 


2 


45 


3 


46 


4 


47 


0 


Zero Vector 


+1 


3 


46 


4 


47 


5 


48 


6 


49 


1 


44 


2 


45 


+2 


305 


397 


306 


398 


307 


399 


308 


400 


303 


395 


304 


396 


+3 


125 


139 


126 


141 


127 


143 


128 


145 


123 


135 


124 


137 


+4 


182 


190 


183 


192 


184 


194 


185 


196 


180 


186 


181 


188 


+5 


252 


260 


253 


262 


254 


264 


255 


266 


250 


256 


251 


258 


+6 


490 


498 


491 


500 


492 


502 


493 


504 


488 


494 


489 


496 


7 


590 


598 


591 


600 


592 


602 


593 


604 


588 


594 


589 


596 


8 


678 


686 


679 


688 


680 


690 


681 


692 


676 


682 


677 


684 


-1 


-8 


696 


680 


686 


681 


688 


682 


690 


683 


692 


684 


694 


685 


-7 


590 


596 


591 


598 


592 


600 


593 


602 


594 


604 


595 


606 


-6 


496 


491 


498 


492 


500 


493 


502 


494 


504 


495 


506 


490 


-5 


251 


258 


252 


260 


253 


262 


254 


264 


255 


266 


250 


256 


-4 


188 


182 


190 


183 


192 


184 


194 


185 


196 


180 


186 


181 


-3 


125 


139 


126 


141 


127 


143 


128 


145 


123 


135 


124 


137 


-2 


397 


306 


398 


307 


399 


308 


400 


303 


395 


304 


396 


315 


-1 


4 


47 


5 


48 


6 


49 


1 


44 


2 


45 


3 


46 


0 


Zero Vector 


+1 


45 


3 


46 


4 


47 


5 


48 


6 


49 


1 


44 


2 


+2 
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+3 
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+4 
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184 


194 
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196 
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+5 
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266 
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256 


251 


258 


252 


260 


253 


+6 


492 


502 


493 


504 


488 


494 


489 


496 


490 


498 


491 


500 


+7 


602 


593 


604 


588 


594 


589 


596 


590 


598 
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600 


592 


+8 


681 


692 


676 


682 
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679 
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0 
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694 


685 


696 


680 


686 


681 


688 


682 


690 


683 
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-6 


506 
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143 
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-2 
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3 


46 


4 
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48 
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45 


3 


46 


4 


+2 
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+3 


145 
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125 


139 
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143 


128 


+4 


186 


181 


188 


182 


190 
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504 
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+7 
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592 
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593 
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594 


589 


596 


590 


598 


591 


+8 


690 


681 


692 


676 


682 


677 


684 


678 


686 


679 


688 


680 
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Using the hysteresis comparator outputs flux C 0 and 
torque C T and the stator flux sector S, the proper output 
vector can be chosen to correct the error due to the 
relation in the switching table. The switching 
configuration is made step by step. The selection of a 
voltage vector at each cycle period Te is carried in order 
to maintain the flux and torque within the limits of two 
hysteresis bands. Several switching tables for three-level 
inverter are presented in [10]. 

A new switching table for the 9-level inverter selector 
has been developed in Table III, to achieve an accurate 
control. In order to simplify, the mechanical rotor speed 
will be considered when assigning the voltage vectors 
needed at each one of those sectors. The speed of the 
stator flux linkage vector is given by the modulus of the 
applied voltage vector. Thus, the voltage vectors will be 
chosen according to the rotor speed [11]. Voltage vectors 
with low amplitude will be chosen for lower speeds. 

4. Simulation Results and Discussion 
The proposed scheme is simulated in DTC with 
multilevel control. The induction motor parameters are as 
follows: Rs=4.85Q, Rr=3.805^, Ls=274mH, Lr=274mH, 
Lm=258mH, p=2, J=31g.m2, V=220V, power=1.5kW 
and speed=1420rpm. All simulations have a sample time 
for the control loop of lOOps; the voltage of the DC bus 
is 514V. The requested space voltage vector, demanded 
by the DTC strategy, is assured as shown by the line 
voltage in Figure. 13 Thus, if compared with a seven-level 
DTC strategy [9]; the dv/dt applied to the motor 
terminals is greatly reduced. Furthermore, the 
electromagnetic torque is also improved. It is possible to 
reduce pulsation amplitude by approximately a factor of 
2 when compared to a seven - level DTC strategy. 

To show the effectiveness of the DTC with 7-level and 9- 
level inverters with SVPWM switching technique a 
simulation work has been carried out on induction motor. 





Figure.7: D-Q axis flux response of 9-level Inverter fed 
DTC IM drive. 




Figure. 8: Phase current response of 7-level Inverter fed 
DTC IM drive . 




Figure. 6: D-Q axis flux response of 7-level Inverter fed 
DTC IM drive 



Figure. 9: Phase current response of 9-level Inverter fed 
DTC IM drive 
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Figure. 10: Torque error for 7-level inverter fed DTC IM 
drive 




Figure. 1 1: Torque error for 9-level inverter fed DTC IM 
drive 





Figure. 14: Stator line voltage harmonic spectrum of 9- 
level inverter 




Figure. 15: Stator line voltage harmonic spectrum of 9- 
level inverter 




Time(sec) 



Figure. 16: Stator flux error of 7 -level Inverter fed 
DTC IM drive 



Figure. 12: Stator line voltage of 7-level inverter 




Time(sec) 



Figure. 13: Stator line voltage of 9-level inverter 




Figure. 17: Stator flux error of 9-level Inverter fed 
DTC IM drive 
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In figure 6 and 7 D-Q axis flux response of the 7 and 9 - 
level inverters are compared. It is seen that the 
performance of the 9-level inverter fed DTC IM drive has 
lower ripple, so the proposed system is superior to 
control the flux with reduced ripple content. Figure 8 and 
9 illustrate the phase current performance of 7-level and 
proposed inverter system, ,the phase current response 9- 
level inverter is fast recovery is obtained .Figure 10 and 
1 1 shows the torque errors of the two systems, from the 
simulation results 7-level inverter torque error is 0.5 N-M 
and 9-level inverter fed DTC IM drive torque error is 
reduced to 0.09 N-M, so from the proposed system 
torque ripple minimization is possible. The stator line to 
line voltages of 7-level and 9-level inverter system are 
illustrated in Figure 12 and 13. 

Figures 14, 15 shows the Stator line voltage harmonic 
spectrum of 7 and 9-level inverter system. The 7-level 
inverter THD is 7.81% and proposed system THD is 
reduced to the 5.99%. From the analysis 9-level inverter 
has reduced THD; it indicates that dv/dt applied to the 
motor terminals is greatly reduced. Furthermore, the 
electromagnetic torque is also improved. Figure 16 and 
17 shows the stator flux errors response of the two 
systems (7 and 9-level inverters). 

From the above discussion, the proposed DTC IM drive 
system behavior is good, even in extreme conditions like 
the reverse speed reference with nominal load torque 
applied. Reduction in ripple is observed in both 
electromagnetic torque and flux is due to the use of 
hysteresis controllers. 

5. Conclusions 

A nine-level inverter system for DTC induction motor 
drive is presented. DTC strategy of a nine-level Multi 
Point Clamped (MPC) voltage source inverter fed 
induction motor is implemented using Matlab/Simulink . 
Nine-level inverter has 729 switching states and there are 
106 effective vectors are possible. The proposed DTC IM 
drive scheme is capable for enough degrees of freedom to 
control both electromagnetic torque and stator flux with 
very low ripple. 

The following points are most attractive features of 
multilevel inverter fed DTC IM drive. 

• The output voltages with extremely low 

distortion and lower dv/dt . 

• They can operate with a lower switching 
frequency. 

• As the number of levels increased 

the %THD in the motor phase voltage 
decreased,. 

• The number of level increased the torque 
ripple is reduced to minimum and the stator 
flux ripple is also minimized. 

• For the proposed system has stator flux 
trajectory response is a circle and answer 
the response is faster as compared to the 7- 
level inverter. 




Nomenclature 

<j> s — Stator flux linkage 

U s = Stator voltage expressed in the stator reference 
frame 

R s = Stator resistance 
I s = Stator current 
At = Standard deviation 

T e = Instantaneous value of the electromagnetic torque 

p = Derivation operator 

ip s - Imaginary component of stator current 

i as = Real component of stator current 

® ^ = D-axis stator flux linkages 

® = Q- axis stator flux linkages 

T ref = Torque reference 

X P sre j = stator flux reference 
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Abstract 

This paper presents a road map for the design and 
development of sustainable hybrid High Brightness Light 
Emitting Diode (HB LED) illumination systems with 
maximum power point tracking (MPPT). The proposed 
hybrid system design presents the foundation for future 
sustainable high efficiency illumination systems. The 
system design presented in this paper individually 
controls an array of HB LEDs. Each HB LED operates 
in a different mode defined by the user. Photovoltaic 
(PV) panels-based systems are used as the primary source 
of energy in this system while the electric grid is used as 
the backup source that supplies power to the HB LEDs 
only when the power of the PV system is not sufficient to 
supply all the system HB LEDs and the PV system 
batteries are depleted below the minimum level set in the 
system design. A matlab/ Simulink model is proposed for 
studying the system behaviors and operation modes. 

Keywords: MPPT, High Brightness LED; Hybrid; 
Matlab/ Simulink; DC Chopper 

1. Introduction 

The demands for utilizing alternative power sources have 
increased due to rising oil prices and more stringent 
environmental regulations. Alternative energy and its 
applications have been heavily studied for the last decade, 
and solar energy is the preferred choice in many 
applications [1]. Among solar energy applications, the 
photovoltaic (PV) technology has received much 
attention, and it is being used in many applications [2&3]. 
This paper presents a sustainable hybrid high brightness 
LED Illumination System that can be used to replace 
current illumination systems in order to improve the 
system efficiency and reliability. In the proposed system, 
a Single-ended primary-inductor converter (SEPIC) DC- 
DC converter is used to deliver solar energy via PV cell 
modules to a battery bank in charging mode during the 
daytime. At night, it drives an LED lighting system. 

The main reason for choosing HB LEDs as the 
illumination source in this system is due to their high 
efficiency as compared with incandescent and fluorescent 
lamps. A second factor is that the world of DC 



Application is fast expanding with the requirements of 
data storage and high efficiency requirements in supply 
distributions. Compact fluorescent lamps are becoming 
obsolete, and replaced with HB LED lamps driving the 
efficiency of lighting systems. The system presented in 
this paper expands on the continued growth of DC power 
distribution in buildings, and driving this trend toward that 
end LED-based lighting has been cited as a factor. 

2. Hybrid LED Illumination 

The principal motivation for this hybrid system is the 
clear shift to DC systems with more use of alternative 
energy sources. The AC vs. DC battle raged when 
Edison promoted DC power while George Westinghouse 
felt that AC was the way to go. As we know, AC won the 
battle, since it was so much easier to step the voltage up 
and down using transformers, and higher voltages greatly 
reduce resistive loss. This paper explores the continued 
growth of DC power, and the force driving the trend; 
LED-based lighting has been cited as one major driver. 
Since photovoltaic panels are used to power the 
illumination system, there is a need for a second source 
of power for the system when the sun is out for a period 
of time beyond the capability of the batteries to serve as a 
backup source. Some products that are appearing in the 
market use HB LED with build-in converter that can be 
connected directly to the AC line. These products take 
advantage of the high efficiency of HB LEDs, but they 
do not use a sustainable source as the main supply for 
those LEDs. The authors present this new hybrid HB 
LED-based illumination system design with automatic 
transfer switch as given in figure 1. This hybrid system 
uses solar as the primary source of energy, and it 
switches to the AC line only when the primary source 
can’t supply the required power to the system. 

In the system given in figure 1, LEDs are powered from 
the solar panel during the daytime as the primary source, 
with the battery in charge mode, and they operate from 
the battery at night. If the battery is fully discharged, the 
HB LED System operates from the AC line as the 
secondary source. The design of DC-DC converters, DC- 
AC rectifiers, and DC- AC inverters are well understood 
and can be followed using senior level or graduate power 
electronics textbooks [4-8] 
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is operated by means of a controller circuit, i.e. an FPGA. 
Figure 5 presents the basic control system layout used for 
this DC illumination control strategy that gives the user 
more flexibility and control than the flexibility and 
control level available for standard AC powered 
illumination used in incandescent or fluorescent lamps. 



U 
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Figure 1 : Hybrid High Brightness LED Illumination System 



Figure 3: AC and DC Contactors Design to Control the LED 
Illumination System Switching 



Since LEDs can handle voltage fluctuations, the system 
can be simplified by using a full wave bridge rectifier 
instead of the AC-DC converter. The simplified system 
is given in figure 2 where components can be classified 
into two different categories. The first one is high 
tension parts which consist of photo cells and high power 
LED, and the other is the low tension side which consists 
of L, C, and battery. The operation of this buck-boost 
converter has two operation modes: the first is the buck 
mode for day time battery charging operation, and the 
second is the boost mode for evening time lighting usage. 




Figure 2: Simplified Hybrid Illumination 



This simplified system is sustainable for many 
applications emerging due to the continued growth of DC 
power distribution in buildings, and the principal force 
behind this trend is LED-based lighting which we have 
cited as one major driver. The main challenges in 
building such a system include analysis of a switching 
circuit that controls alternating of power between the AC 
and DC distribution systems, being fed from Solar panel 
and switching between the two sources. Figure 3 
presents one solution for this problem where an AC 
contactor and a DC contactor with mechanical interlock 
are used, and figure 4 presents the relay control for the 
mechanically interlocked contactors. 

The illumination system is designed to be used to 
individually control an array of HB LEDs. Each HB LED 
operates in a different mode that is independently defined 
by the user. The defined sequence of LED illumination 




PLA Relay 
Pick-up Voltage 
= 40.5V 



Figure 4: Relay Design used to Control the Hybrid Illumination 
Switching 




Figure 5: FPGA-Controlled Illumination 

The FPGA is programmed using Hardware Description 
Language (HDL). In this case VHDL is used. [VHDL is 
a hardware description language used in electronic design 
automation to describe digital and mixed systems and 
integrated circuits.] The pattern of the LEDs defined by 
the user can be modified using this hardware language. 
Basys Spartan 2 FPGA is used for this application due to 
its flexibility, low cost, and ease of use. 

The application software used in our case is Adept which 
is a powerful program that allows configuration and data 
transfer with Xilinx logic devices, and is used as an 
interface between Xilinx and the Spartan 2 FPGA. For 
maximum intensity from the LEDs, the typical forward 
voltage of 3.9V, with forward current 700mA, is supplied 
to the FPGA. A personal computer provided with Xilinx 
software environment is used for programming the 
sequence of the LEDs defined by the user. The 
programming is done again using VHDL. The developed 
software was tested using Digilent FPGA. 
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High brightness LEDs [9&10] can be driven at currents 
from hundreds of mA to more than an ampere, compared 
with the tens of mA for other LEDs; however few of the 
HB LEDs can produce over a thousand lumens. Since 
overheating is destructive, the HP LEDs may need to be 
mounted on a heat sink to allow for heat dissipation. If 
the heat from an HB LED is not removed using a heat 
sink, the device will burn out in seconds. 

A single HB LED can often replace an incandescent bulb 
in a torch, or be set in an array to form a powerful HB 
LED system. LEDs can operate on AC power without the 
need for a DC converter. Each half-cycle part of the LED 
emits light with the other half cycle being dark, and this 
is reversed during the next half cycle. The efficacy of this 
type of HB LED is typically around 40 LM/W. A large 
number of LED elements in series may be able to operate 
directly from line voltage. In 2009 Seoul Semiconductor 
released a high DC voltage capable of being driven from 
AC power with a simple controlling circuit. The low 
power dissipation of these LEDs gives them more 
flexibility than the original AC LED design. In this 
project the HB LEDs are powered from a DC source. 

To control the LED from the FPGA, a driver circuit is 
needed to boost the power of the FPGA output. An LED 
driver circuit is an electric circuit used to power a light- 
emitting diode or LED. The circuit consists of a voltage 
source powering a current limiting resistor and an LED 
connected in series. The HB LEDs used in our design 
have a constant current of 700mA and a supply voltage of 
+3V. As the current has to be amplified to 700mA for 
each LED, two transistors are connected together so that 
the current amplified by the first is amplified further by 
the second transistor. The overall current gain is equal to 
the two individual gains multiplied together, i.e h FE = 
hpEi x h FE2 , where h Fm and h FE2 are the gains of the 
individual transistors. 

3. Sustainable PV-Powered Illumination 

Since photovoltaic panels are used to power the HB LED 
illumination systems, there is a need for the development 
of the DC converter system to power the LEDs. To 
satisfy this requirement, a forward converter with PV- 
based LED applied in lighting systems was used. In the 
proposed DC supply system, the SEPIC is used to deliver 
solar energy via PV cell modules to a battery bank in 
charging mode during the daytime. During the nighttime, 
the converter drives an LED lighting system. 

Figure 7 illustrates the principle PV-SEPIC- LED circuit 
applied in street illumination, where the produced PV 
voltage is stored in a battery bank throughout the 
charging unit during the daytime, at night the discharger 
activates and the LEDs are energized with appropriate 
voltage throughout the step-up transformer and full 
bridge rectifier. 

SEPIC circuits find widespread application when the 
input voltage fluctuates above and below average value, 
while the output voltage must be kept at a constant value 
with minimum tolerances. One of most important 
applications of SEPIC circuit is the integration with the 
photovoltaic system (PV system) and illumination load of 
series and parallel connected LEDs. 

The SEPIC converter is a DC/DC converter topology that 
provides a positive regulated output voltage from an 



input voltage that varies from above to below the output 
voltage. This type of topology is needed when the 
voltage from an unregulated input power source such as 
solar where the sun irradiation, temperature and weather 
changes directly affect the generated output voltage. The 
standard SEPIC topology [11&12] requires two inductors 
in additional to step-up transformer, making the power- 
supply footprint quite large. 

Photovoltaic (PV) cells are used to convert the sunlight 
into electrical energy. On the other hand, it is also an 
important issue to save the energy demand and increase 
the energy efficiency [13-15]. High brightness light 
emitting diodes (LEDs) [16&17] are becoming more 
widespread for the lighting applications such as 
automobile safety and signal lights, traffic signals, street 
lighting and so on. LED power circuitry is discussed 
thoroughly in the literature [18]. 

In lighting applications with solar energy, the charger is 
adapted to convert the solar irradiations for storing in the 
battery during the daytime. In the nighttime, a discharger 
is used to release energy in the battery and drive the LED 
lighting system. Low-power DC-DC converters can be 
used for the charger and discharger mode. Since the PV 
voltage from the solar panel is unstable, the buck-boost 
converter is more suitable for the charger circuit. This 
converter can also be used in the discharger circuit. 

4 . DC-DC Converter MPPT Control Design 
and Analysis Using Computer Simulation 

Computer simulation is an important tool for future 
illumination systems design. The HB LED-based 
illumination system was simulated using a 

matlab/Simulink environment [19]. In this section, we 
present the complete simulation model including the 
solar PV module, Maximum Power Tracking module 
(MPPT), DC-DC boost converter and storage unit 
(Battery ). These modules are going to be described as 
follows: 



4.1 Characteristics of PV Array 

Basically, PV cell is a P-N semiconductor junction that 
directly converts light energy into electricity. It has the 
equivalent circuit shown in Figure 1 as represented in [20- 
22 ]. 

Rs Io 





| Vo 



■o 



Figure 6: Equivalent circuit for PV cell 



Where I ph represents the cell photocurrent; R p and R s are 
the intrinsic shunt and series resistance of the cell 
respectively; Id is the diode saturation current; V 0 and I 0 
are the cell output voltage and current respectively. 

The following are the simplified equations describing the 
cell output voltage and current: 



Vo = 



A.K.Tc 



In 



Iph + Id — Io 

Io 



— Rs.Io 



( 1 ) 
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( q.Vo/Ns ^ 
* A.K.Tc _ | 



l-V performance 



Id = Io 



^ 3 q-Eg f 1 1_ 

B . K v Tr Tc 

• V 



.Try 
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Where, K- Boltzmann constant; N p and Ns are the 
number of parallel and series connected cells respectively; 
E g is the band gap of the semiconductor; T c and T r are the 
cell and the reference temperature respectively in Kelvin, 
A and B are the diode ideality factors with values varies 
between 1 and 2; <D n is the normalized insulation; I sc is 
the short circuit current given at standard condition; I t 
and I or are constants given at standard conditions. 

4.1 Photovoltaic I-V Performance 

In order to study the I-Y performance of the PV circuit 
and to look for appropriate dc chopper for boosting up 
the output voltage to predetermined value it is necessary 
to illustrate the obtained PV voltage and current for 
boost chopper according to specifications given in table 1 
at reference irradiation 1000W/m2. 

The PV Array voltage can be obtained by multiplying the 
module voltage and current 
represents number of series 
modules respectively. 



by N sm and N pm that 
and parallel connected 



Table 1: Data specification for PV Array. 
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a) Matlab/Simulink module for PV Array 
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b) I-V Performance of PV module. 



Power performance 




c) Power performance of proposed module 
Figure 7: PV Array module & main charaterestics. 

Figure 7 illustrates the proposed PV array with R load 
presenting the equivalent resistance of complete LED set, 
where the obtained results for different variation levels 
are presented. From these performances it is shown that 
the total output PV voltage and current varies according 
to irradiation level with approximated 65W maximum 
power at G=1000W/m2. 

4.2 Maximum power point derivation 

In order to operate the module at maximum extracted 
power, it is necessary to calculate the coordinates of the 
maximum power point (V M pp, Impp)- For this, and to 
simplify the model in Simulink, the coordinates of the 
maximum power point are given by the following 
equations: 

Impp = Iph - Id.(e Vmpp/VT - 1) (5) 

Iph — Impp / Voc/Vt 



Vmpp = VT.ln 1 + - 



Iph 



( 6 ) 



Where Vt and Voc are the thermal voltage 
and the cell open circuit voltage respectively 
and given by: 

K.Tc 



Vt = 



V oc = V t . In 



(V) 



— +1 

Id 



The maximum power that can be obtained can be 
expressed as follows: 



MPP 



= V 



MPP 



.1 



MPP 



(9) 
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Obtaining ng the maximum obtained power form the 
PV for given solar irradiation requires controlling the duty 
cycle D of DC Buck-Boost Chopper toward operating the 
chopper in mode where the extracted from the system 
power will be maximum. There are several approaches 
used for realizing MPPT operation mode, one of the most 
popular and simple approach de is called Perturbation & 
Observation model (P&O MPPT model). 

The mathematical equations realized MPPT approach are 
expressed as follows: 

Vo = Vpv. — - — 

1 - D 

Po = Vo 2 /RI 

dPo = Po(k + 1) - Po(k); 

dVo = Vo(k + 1) - Vo(k); (io) 

(dPo/dVo) = ? -> D = ? 

[> 0 — » D = D - AD 
[< 0 — » D = D + AD 
[>0— »D = D + AD 
< 0 — > D = D - AD 



dPo 



> 0 -> dVo 



< 0 -> dVo 



While figure 8 illustrates the Power performance and the 
flowchart procedure for obtaining the appropriate duty 
cycle needed to operate the chopper at maximum obtained 
power. 




The proposed mathematical equations and corresponds 
flowchart can be summarized with operation logic as 
follow: 

• (dPo/dVo)<0 ^ operation at the right hand side of 
the curve (fig. 8 a), for D-> Dmin , then V must be 
decreased ( moving toward decrease). 

• (dPo/dVo>0) 4 operation at the left hand side of 
the same curve , for D-> Dmax^l). Then V must be 
increased (moving toward increase). 

4.3: Boost DC Chopper 

The output PV voltage in most of applications is requred 
to be boost up by applying Boost DC chopper [6] as 
shown in figure 9, where the principle electrical circuit 
and the operation mode at continuous current mode 
(CCM) are illustrated. 





b) Flowchart for tracking procedure 
Figure 8: P&O MPPT tracking method. 
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c) Mode#2: IGBT=OFF 




d) Load voltage & current at two operation modes 
Figure 9: Equivalent circuit for boost converter operating in CCM 
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4.4: Matlab Simulation of Proposed Model 

The complete electrical circuit for LED illumination 
system is simulated in mat lab/simulink environment and 
illustrated in 10 




D1 





c) MPPT module 

Fig. 10: Complete Simulation model of PV- LED illumination system. 



4.5: Simulation Results 

There are various PV and chopper paramtetrs that can he 
described; some of these are illustrated in figure 11 as 
follows: 

• Solar irradiation : Figure 11(a) illustrates the time 
profile of solar irradiation levels for various periods 
of daytime in order to show the behaviors of the 
circuit according to these levels. 

• Chopper duty cycle : Figure 11(b) illustrates the time 
profile of calculated dudty cycle , where its shown 
that for each irradiation level the duyt cycle increases 
which cause further increase until maximum obtaied 
power is achieved. 

• PV array current'. Figure 11(c) illustrates the time 
profile of PV current, where its shown that high 
irradiation produces large current and power . 



Irradiation vs simulation time 




Simulation interval 



a) Solar irradiations at various daytime intervals. 



Duty cCycle vs simulation time 




b) Chopper duty cycle. 



PV Array current vs simulation interval 




c) Array current of The PV pannel 




d) PV voltage, current and power 
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Array & output voltage 




e) PV and chopper voltage 
Fig. 1 1 : Circuit performances at various irradiations 



• PV voltage , current and power. Figure 11(d) 
illustrates the time profile of PV voltage , current and 
power , where its shown how the power being 
affected by increasing the irradiation level. 

• Input-output voltage : Figure 11(e) illustrates the 
time profile of PV input voltage and output chopper 
voltage, where its shown how the input voltage varies, 
while the output voltage is boosted up and fixed at 
appropriate for the elimination system value. 

5. Conclusion 

A Sustainable hybrid MPPT-Controlled HB LED-Based 
Illumination System was developed. This illumination 
system can be used for many lighting applications since it 
is more efficient and reliable than existing traditional 
lighting systems based on incandescent or fluorescent 
lamps. While the use of HB LEDs adds to the initial cost 
of the system, it will be paid off in the long term due to 
their higher reliability, flexible control, and long life time. 
A simulink model for this illumination system was built 
in order to study the system behaviors at various 
irradiation level, the effect of MPPT system and 
generated duty cycle on chopper operation, which in turn 
causes better utilization, effecienecy enehancement and 
loss reduction. 
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Abstract 

A sophisticated space vector modulation scheme for a 
nine-level inverter system for dual-fed induction motor 
drive, using only the instantaneous sampled reference 
signals is presented in this paper. The dual-fed structure 
is realized by opening the neutral-point of the 
conventional squirrel cage induction motor. The nine- 
level inversion is obtained by feeding the dual-fed 
induction motor with symmetrical four-level inverter 
from one end and symmetrical three-level inverter from 
other end. The proposed space vector pulse width 
modulation technique does not require the sector 
information and look-up tables to select the appropriate 
switching vectors. The inverter leg switching times are 
directly obtained from the instantaneous sampled 
reference signal amplitudes and centers the switching 
times for the middle space vectors in a sampling time 
interval, as in the case of conventional space vector pulse 
width modulation. 



Keywords: dual-fed induction motor, middle space 
vectors, nine-level inverter, sampled sinusoidal reference 
signals, space vector PWM. 



Nomenclature 

Edc 



Ts 

Ea3A5> E B 3B5? E C 3C5 
Emax, Emin 



Eoffsetl, Eoffset2 

E an? E bn , E cn 



ti, t 2 



toi? t()2 

Tx-cross(x=a, b, c) 



DC link voltage of an equivalent 
conventional single two-level 
Inverter 
Sampling time 
Motor phase voltages 
The maximum and minimum 
magnitudes of the three sampled 
sinusoidal reference signals 
common offset voltages 
The modified reference signals 
after addition of Eoffsetl 
The time durations for the middle 
inverter space vectors 
The time durations for the start and 
end inverter space vectors 
The time interval, at which the x- 
phase reference signal, E xn 
crosses the triangular carrier signal 



T*as, T*bs, T*cs 

Tfirst_cross, 

Tsecond_cross, 

Tthird_cross 



T 0 

Tmiddle 



The time equivalents of the voltage 
magnitudes 



The time intervals, at which the 
sinusoidal reference signals cross 
the triangular carrier signals for the 
first time, second time and third 
time 

The time interval of the start and 
end space vector 

The time interval for the middle 
inverter space vectors 



1. Introduction 

The two most widely used pulse width modulation 
(PWM) schemes for multilevel inverters are the carrier- 
based sine-triangle PWM (SPWM) scheme and the space 
vector PWM (SVPWM) scheme. The SPWM schemes are 
more flexible and simpler to implement, but the maximum 
peak of the fundamental component in the output voltage 
is limited to 50% of the DC link voltage [1, 2]. The 
maximum peak of the fundamental component in the 
output voltage obtained with space vector modulation is 
15% greater than with the sine-triangle modulation 
scheme [2-5]. But the conventional SVPWM scheme 
requires sector identification and look-up tables to 
determine the timings for various switching vectors of the 
inverter, in all the sectors [3, 4]. This makes the 
implementation of the SVPWM scheme quite 
complicated. It has been shown that, for two-level 
inverters, a SVPWM like performance can be obtained 
with a SPWM scheme by adding a common mode voltage 
of suitable magnitude, to the sinusoidal reference signals 
[4-6]. 

The SPWM scheme, when applied to multilevel inverters, 
uses a number of level-shifted carrier signals to compare 
with the sinusoidal reference signals [7]. The SVPWM for 
multilevel inverters [8 - 11] involves mapping of the outer 
sectors to an inner sub-hexagon sector, will be very 
complex, as a large number of sectors and inverter vectors 
are involved. 
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A modulation scheme is presented in [12], where a fixed 
common mode voltage is added to the reference signal 
throughout the modulation range. It has been shown [13] 
that this common mode voltage addition will not result in 
a SVPWM-like performance, as it will not centre the 
middle space vectors in a sampling interval. The common 
mode voltage to be added in the reference phase voltages, 
to achieve SVPWM-like performance, is a function of the 
modulation index for multilevel inverters [13]. A 
SVPWM scheme based on the above principle has been 
presented [14], where the switching time for the inverter 
legs is directly determined from sampled sinusoidal 
reference signal amplitudes. But it involves region 
identifications based on modulation indices. While this 
SVPWM scheme works well for a three-level PWM 
generation, it cannot be extended to multilevel inverters of 
levels higher than three, as the region identification 
becomes more complicated. A carrier-based PWM 
scheme has been presented [15], where sinusoidal 
references are added with a proper offset voltage before 
being compared with carriers, to achieve the performance 
of a SVPWM. The offset voltage computation is based on 
a modulus function depending on the DC link voltage, 
number of levels and the sinusoidal reference signal 
amplitudes. A SVPWM scheme is presented [18], where 
the switching time for the inverter legs is directly 
determined from sampled sinusoidal reference signal 
amplitudes for five-level inverter where two three-level 
inverters feed the dual-fed induction motor. A carrier 
based SPWM scheme is presented [19, 20] for five-level 
and nine-level inverters. Qualitative SVPWM schemes are 
presented for five-level inverters in [21, 22]. 

The objective of this paper is to present an 
implementation scheme for PWM signal generation for 
nine-level inverter system for dual-fed induction motor, 
similar to the SVPWM scheme. In the proposed scheme, 
the dual-fed induction motor is fed with symmetrical 
four-level inverter from one end and symmetrical three- 
level inverter from other end. The PWM switching times 
for the inverter legs are directly derived from the sampled 
amplitudes of the sinusoidal reference signals. A simple 
way of adding an offset voltage to the sinusoidal 
reference signals, to generate the SVPWM pattern, from 
only the sampled amplitudes of sinusoidal reference 
signals, is explained. The proposed SVPWM signal 
generation does not involve checks for region 
identification, as in the SVPWM scheme presented in 
[14]. Also, the algorithm does not require either sector 
identification or look-up tables for switching vector 
determination as are required in the conventional 
multilevel SVPWM schemes [10, 11]. Thus the scheme is 
computationally efficient when compared to conventional 
multilevel SVPWM schemes, making it superior for real- 
time implementation. 

2. Nine-level inverter scheme for the dual- 
fed induction motor 

The power circuit of the proposed drive is shown in 
Figure 1. A symmetrical four-level inverter, Inverter-A 
and a symmetrical three-level inverter, Inverter-B feed 
the dual-fed induction motor. The inverter-A is composed 
of three conventional two-level inverters INV-1, INV-2 



and INV-3 in cascade. The Inverter-B is composed of 
two conventional two-level inverters INV-4 and INV-5 in 
cascade. The DC link voltages of INV-1, INV-2, INV-3, 
INV-4 and INV-5 are (2/8)Edc, (2/8)Edc, (2/8)Edc, 
(l/8)Edc and (l/8)Edc respectively, where Edc is the DC 
link voltage of an equivalent conventional single two- 
level inverter drive. 

The leg voltage E A 3 n of phase-A attains a voltage of 
(2/8)Edc if (i)The top switch S 3i of INV-3 is turned on 
(Figure 1) and (ii) The bottom switch S 2 4 of INV-2 is 
turned on. The leg voltage E A 3 n of phase-A attains a 
voltage of (4/8)Edc if (i) the top switch S 3i of INV-3 is 
turned on (ii) The top switch S 2 i of INV-2 is turned on 
and (iii)The bottom switch Si 4 of INV-1 is turned on. The 
leg voltage E A3n of phase-A attains a voltage of (6/8)Edc 
if (i) the top switch S 3i of INV-3 is turned on (ii) The top 
switch S 2 i of INV-2 is turned on and (iii)The top switch 
Sn of INV-1 is turned on. The leg voltage E A3n of phase- 
A attains a voltage of zero volts if the bottom switch S 34 
of the INV-3 is turned on. Thus the leg voltage E A3n 
attains four voltages of 0, (2/8)Edc, (4/8)Edc and 
(6/8)Edc, which is basic characteristic of a 4-level 
inverter. Similarly the leg voltages E B 3 n and E C 3 n of 
phase-B and phase-C attain the four voltages of 0, 
(2/8)Edc, (4/8)Edc and (6/8) Edc. 

The leg voltage E A5n ’ of phase-A attains a voltage of 
(l/8)Edc if (i)The top switch S 5 i of INV-5 is turned on 
and (ii) The bottom switch S 44 of INV-4 is turned on. The 
leg voltage E A5n ’ of phase-A attains a voltage of (2/8)Edc 
if (i) The top switch S 5 i of INV-5 is turned on and (ii) 
The top switch S 4) of INV-4 is turned on. The leg voltage 
E A5n ’ of phase-A attains a voltage of zero volts if the 
bottom switch S 54 of the INV-5 is turned on. Thus the leg 
voltage E A5n ’ attains three voltages of 0, (l/8)Edc and 
(2/8)Edc, which is basic characteristic of a 3 -level 
inverter. Similarly the leg voltages E B5n ’ and E C 5 n ’ of 
phase-B and phase-C attain the three voltages of 0, 
(l/8)Edc and (2/8)Edc. 

Thus, one end of dual-fed induction motor may be 
connected to a DC link voltage of either zero or (2/8)Edc 
or (4/8)Edc or (6/8)Edc and other end may be connected 
to a DC link voltage of either zero or (l/8)Edc or 
(2/8)Edc. When both the inverters, Inverter-A and 
Inverter-B drive the induction motor from both ends, nine 
different levels are attained by each phase of the 
induction motor. The nine levels generated for phase-A 
are shown in Table 1. 



TABLE 1 

THE NINE LEVELS REALIZED IN THE PHASE-A WINDING 



Leg- 
voltage 
of phase A, 

E A3n 


Leg- 
voltage 
of phase A, 

E A 5n' 


Motor phase voltage 

E A3A 5 = E A3n - E A5n - 


Level 


0 


(2/8) Edc 


-(2/8) Edc 


Level 1 


0 


(l/8)Edc 


-(l/8)Edc 


Level 2 


0 


0 


0 


Level 3 


(2/8) Edc 


(1/8) Edc 


(l/8)Edc 


Level 4 


(2/8) Edc 


0 


(2/8)Edc 


Level 5 


(4/8) Edc 


(1/8) Edc 


(3/8)Edc 


Level 6 


(4/8) Edc 


0 


(4/8) Edc 


Level 7 


(6/8) Edc 


(1/8) Edc 


(5/8) Edc 


Level 8 


(6/8) Edc 


0 


(6/8) Edc 


Level 9 
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INVERTER A 




Figure 1 . Schematic circuit diagram of the proposed 9-level inverter drive scheme. 



3. Proposed SVPWM in linear modulation 
range 

For two-level inverters, in the SPWM scheme, each 
sinusoidal reference signal is compared with the triangular 
carrier signal and the individual phase voltages are 
generated [1]. To attain the maximum possible peak 
amplitude of the fundamental phase voltage, a common 
offset voltage, Eoffsetl is added to the sinusoidal 
reference signals [5, 12], where the magnitude of Eoffsetl 
is given by 

Eoffsetl = -(E max+ E min) / 2 (1) 

Where Emax and Emin are the maximum and minimum 
magnitudes of the three sampled sinusoidal reference 
signals respectively, in a sampling time interval. The 
addition of Eoffsetl, results in the active space vectors 
being centered, making the SPWM scheme equivalent to 
the SVPWM scheme [3]. In a sampling time interval, the 
sinusoidal reference signal which has lowest magnitude 
crosses the triangular carrier signal first, and causes the 
first transition in the inverter switching state and the 
sinusoidal reference signal, which has the maximum 
magnitude, crosses the triangular carrier signal last and 
causes the last switching transition in the inverter 
switching states in a two-level SVPWM scheme [5, 13]. 
For nine-level inverter, the modified sinusoidal reference 
signals (E AN ^ E* B n and E C n) after addition of offset 



voltage Eoffsetl, are shown in figure 2 along with eight 
triangular carrier signals T1 to T8. Each time a sinusoidal 
reference signal crosses the triangular carrier signal, it 
causes a change in the inverter switching state. The 
changes in phase voltage and their time intervals are 
shown in Figure 3 in a sampling time interval Ts. The 
sampling time interval Ts can be split into four time 
intervals t 0 i, ti, t 2 and t 02 . The time intervals t 0 i and to 2 are 
the time durations for the start and end inverter space 
vectors respectively and the time intervals ti and t 2 are the 
time durations for the middle inverter space vectors 
(active space vectors), in a sampling time interval Ts. It 
should be observed from Figure 3 that the middle space 
vectors are not centered in a sampling time interval Ts. 
Because of the level- shifted eight triangular carrier signals 
(Figure 2), the first crossing and the last crossing of the 
sinusoidal reference signal cannot always be the minimum 
and the maximum magnitude of the three sampled 
sinusoidal reference signals, in a sampling time interval. 
Thus the offset voltage, Eoffsetl is not sufficient to center 
the middle inverter space vectors, in a multilevel PWM 
system (Figure 3). Hence an additional offset (offset2) has 
to be added to the sinusoidal reference signals of Figure 2, 
so that the middle inverter space vectors can be centered 
in a sampling time interval, same as a two-level SVPWM 
system [3]. 
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4. Determination of the offset voltage for a 
nine-level inverter 

Figure 2 shows modified sinusoidal reference signals and 
eight triangular carrier signals used for PWM generation 
for nine-level inverter. The modified sinusoidal reference 
signals are given by 

e *an = e an + Eoffsetl 

e *bn ~ e bn + Eoffsetl (2) 

e *cn = e cn + Eoffsetl 

where E AN , E B n and E C n are the sampled amplitudes of 
sinusoidal reference signals during the current sampling 
time interval and Eoffsetl is calculated from equation (1). 
The time interval, at which the A-phase sinusoidal 
reference signal, E an crosses the triangular carrier signal, 
is termed as Ta-cross (Figure 4). Similarly, the time 
intervals, when the B -phase and C-phase sinusoidal 
reference signals, E* BN and E* C n cross the triangular 
carrier signals, are termed as Tb-cross and Tc-cross 
respectively. Figure 4 shows a sampling time interval 
when the A-phase sinusoidal reference signal is in the 
triangular carrier region T7 while the B -phase sinusoidal 
reference signal and C-phase sinusoidal reference signal 
are in carrier region T8 and T1 respectively. As shown in 
Figure 4, the time interval, Ta-cross, at which the A- 
phase sinusoidal reference signal crosses the triangular 
carrier signal is directly proportional to the phase voltage 
amplitude, (E* AN - 3Edc/8). The time interval, Tb-cross, 
at which the B -phase sinusoidal reference signal crosses 
the triangular carrier signal, is proportional to (E BN + 
4Edc /8) and the time interval, Tc-cross, at which the C- 
phase sinusoidal reference signal crosses the triangular 
carrier signal, is proportional to (E CN ). Therefore 

Ta - cross = ( E -3Edc/8)\ — — — | = T * as - 3 Ts 

\Edd 8) 

Tb - cross = (E* bn +4 Edc / 8) — } = T*bs + 4Ts 

\Edd 8 j 

Tc - cross = ( E * CN ^ j g J ~ ^ * cs ^ 

Where T*as, T*bs and T*cs are the time equivalents of the 
voltage magnitudes. The proportionality between the time 
equivalents and corresponding voltage magnitudes is 
defined as follows [6]: 

( Edc/8)/Ts = E* AN IT* as 
(Edc 18) ITs = E * BN IT *bs 
(Edc / 8) / Ts = E* cn IT*cs 

(Edc IS)/Ts = Eoffsetl / Toffsetl (4) 

The time interval, at which the sinusoidal reference 
signals cross the triangular carrier signals for the first 
time, is termed as Tfirst_cross. Similarly, the time 
intervals, at which the sinusoidal reference signals cross 
the triangular carrier signals for the second and third time, 
are termed as, Tsecond_cross and Tthird_cross 
respectively, in a sampling time interval Ts. 

Tfirst_ cross = min(T a - cross , Tb - cross , Tc- cross) 

T second _ cross = mid(Ta - cross ; Tb - cross,Tc - cross) 
Tthird_cross= ma x(Ta-cross,Tb-cross,Tc- cross) (5) 



The time intervals, Tfirst_cross, Tsecond_cross and 
Tthird_cross, directly decide the switching times for the 
different inverter voltage vectors, forming a triangular 
sector, during one sampling time interval Ts. The time 
intervals for the start and end space vectors, are t 0 i= 
Tfirst_cross, t 0 2 = (Ts - Tthird_cross), respectively 
(Figure 3). The middle space vectors are 
centered by adding a time offset, Toffset2 to Tfirst_cross, 
Tsecond_cross and Tthird_cross. The time offset, Toffset2 
is determined as follows. The time interval for the middle 
inverter space vectors, Tmiddle, is given by: 

Tmiddle - Tthird _ cross - Tfirst _ cross (6) 

The time interval of the start and end space vector is 

T 0 - Ts- Tmiddle (7) 

Thus the time interval of the start space vector is given by 

T 0 I2 = Tfirst _ cross + Toffsetl 
Therefore 

Toffsetl = T 0 /2- Tfirst _ cross (8) 

In this way, we can obtain offset voltages to be added for 
remaining samples during the time period of sinusoidal 
reference signal and for different modulation indices. 




Figure 2. Modified sinusoidal reference signals and triangular carrier 
signals for a nine-level PWM scheme 
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Figure 3. Inverter switching vectors and their switching time durations 
during sampling time interval Ts 
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Figure 4. Determination of the Ta-cross, Tb-cross and Tc-cross 
during sampling interval Ts 



5. Simulation results and discussion 

The proposed SVPWM scheme is simulated using 
MATLAB environment with open loop E/f control for 
different modulation indices. The speed reference is 
translated to the frequency and voltage commands 
maintaining E/f. The modified three reference sinusoidal 
signals which are added by the total offset voltage to 
make SPWM scheme equivalent to the SVPWM scheme, 
are simultaneously compared with the triangular carrier 
set. A DC link voltage (Edc) of 800 volts is assumed for 



simulation studies. Figure 5(a) shows the total offset 
voltage to be added to sinusoidal reference signals to 
make SPWM equivalent to the SVPWM in the lowest 
speed range which corresponds to two-level mode when 
modulation index is 0.1, figure 5(b) shows the A-phase 
sinusoidal reference signal after offset voltage is added 
and Figure 5(c) shows the motor phase voltage (Easas)- 
Figure 6 to Figure 12 show the motor waveforms in the 
next speed ranges which corresponds to three-level mode 
to nine-level mode when modulation indices are 0.2, 0.3, 
0.42, 0.53, 0.62, 0.74 and 0.85 respectively. It can be 
observed that the motor phase voltage during 9-level 
operation is very smooth and close to the sinusoid with 
lower harmonics. 



6. Conclusion 

A modulation scheme of SVPWM for nine-level inverter 
system for dual-fed induction motor drive is presented. 
The centering of the middle inverter space vectors of the 
SVPWM is accomplished by the addition of an offset 
voltage signal to the sinusoidal reference signals, derived 
from the sampled amplitudes of the sinusoidal reference 
signals. The SVPWM technique, presented in this paper 
does not require any sector identification, as is required in 
conventional SVPWM schemes. The proposed scheme 
eliminates the use of look-up table approach to switch the 
appropriate space vector combination as in conventional 
SVPWM schemes. This reduces the computation time 
required to determine the switching times for inverter 
legs, making the algorithm suitable for real-time 
implementation. 



$ 

0 

> 



V 




“1 — 

V 




— 1— 

V 




r 

V 




V 


1 


VI 


“1 




/\ 


J 


/\ 


L_ 


/\ 


L 


/\ 




/\ 

j 




/\ 





0 0.05 0.1 0.15 0.2 0.25 



Time in Seconds 




100 



Time in Seconds 

M 






0.05 



0.1 0.15 

Time in Seconds 

(0 



0.2 



0.25 



Figure 5. (a)The total offset voltage to be added to sinusoidal reference 
signals (b)The A-phase sinusoidal reference signal after offset voltage 
is added (c) Motor phase voltage when M=0. 1 (2-level operation) 
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Figure 6. (a)The total offset voltage to be added to sinusoidal reference 
signals (b)The A-phase sinusoidal reference signal after offset voltage 
is added (c) Motor phase voltage when M=0.2 (3-level operation) 
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Figure 8. (a)The total offset voltage to be added to sinusoidal reference 
signals (b)The A-phase sinusoidal reference signal after offset voltage 
is added (c) Motor phase voltage when M=0.42 (5 -level operation) 



Volts 


VvWvV 

1 1 1 f 1 1 1 1 


100 

! • 

> 

inn 


'VVWVV- 


-IUU 

1 


1 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.09 0.09 ' I 

Time in Seconds 

(a) 


I 0.01 0.02 0.03 0.04 0.05 0.06 

Time in Seconds 

(a) 


200 

tfl 


1 1 1 [ 1 1 1 1 


400 

200 

00 tuu 


I 1 1 1 1 1 I 


□ o 




O 0 




> 

200 


1 1 1 f 1 1 1 1 


> 200 
400 


: i i i i i : 


[ 


1 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.09 

Time in Seconds 

(b) 


[ 


1 0.01 0.02 0.03 0.04 0.05 

Time in Seconds 

(b) 


200 

| • 
200 


1 1 1 r 1 1 1 1 

1 l l [ l l l l 


200 

! • 
> 200 


'VVy 

: ■ 



0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 

Time in Seconds 

(c) 



0 0.01 0.02 0.03 0.04 0.05 

Time in Seconds 

(0 



Figure 7. (a)The total offset voltage to be added to sinusoidal reference 
signals (b)The A-phase sinusoidal reference signal after offset voltage 
is added (c) Motor phase voltage when M=0.3 (4-level operation) 



Figure 9. (a)The total offset voltage to be added to sinusoidal reference 
signals (b)The A-phase sinusoidal reference signal after offset voltage 
is added (c) Motor phase voltage when M=0.53 (6-level operation) 
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Figure 10. (a)The total offset voltage to be added to sinusoidal reference 
signals (b)The A-phase sinusoidal reference signal after offset voltage 
is added (c) Motor phase voltage when M=0.62 (7-level operation) 



Figure 12. (a)The total offset voltage to be added to sinusoidal reference 
signals (b)The A-phase sinusoidal reference signal after offset voltage 
is added (c) Motor phase voltage when M=0.85 (9-level operation) 
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Figure 1 1 . (a)The total offset voltage to be added to sinusoidal reference 
signals (b)The A-phase sinusoidal reference signal after offset voltage 
is added (c) Motor phase voltage when M=0.74 (8-level operation) 
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Abstract 

This paper presents a simple and novel generalized PWM 
algorithm for n-level multilevel inverter fed direct torque 
controlled induction motor drives. The proposed 
algorithm does not use the angle and sector information 
and hence reduces the complexity involved in the PWM 
algorithm. The proposed algorithm has been developed 
by using the concept of imaginary switching times, which 
are proportional to the instantaneous phase voltages only. 
The proposed algorithm generates the modulating 
waveforms for a n-level inverter, which will be compared 
with the level shifted triangular carriers to generate the 
switching pulses to the inverter. To validate the proposed 
algorithm several numerical simulation studies have been 
carried out and results are presented. From the 
simulation results, it can be observed that the proposed 
algorithm gives reduced harmonic distortion with 
increased number of levels. 

Keywords: Space vector PWM, Multi Level Inverter DTC, 
Induction Motor. 

1. Introduction 

All most the DC drives are replaced by AC drives 
because of its variable speed characteristics. The 
invention of field oriented control (FOC) algorithm has 
been made a renaissance in the high-performance 
variable speed drive applications. The FOC algorithm 
gives the decoupling control of torque and flux of an 
induction motor drive and control the induction motor 
similar to a separately excited dc motor [1]. But, the 
complexity involved in the FOC algorithm is more due to 
reference frame transformations. To reduce the 
complexity involved, a new control strategy called as 
direct torque control (DTC) has been proposed in [2]. A 
detailed comparison is brought between FOC and DTC in 
[3] and concluded that DTC gives fast torque response 
when compared with the FOC. Though, FOC and DTC 
give fast transient and decoupled control, these operate 
the inverter at variable switching frequency due to 
hysteresis controllers. Moreover, the steady state ripples 



in torque, flux and currents are high in DTC. To reduce 
the steady state ripples, discrete space vector modulation 
(DSVM) algorithm has been proposed in [4]. Though, 
DSVM algorithm reduces the torque ripple up to some 
extent, it gives variable switching frequency operation of 
the inverter. 

To reduce the harmonic distortion and to obtain the 
constant switching frequency operation of the inverter, 
nowadays many researchers have focused their interest 
on pulsewidth modulation (PWM) algorithms. A detailed 
survey on various PWM algorithms is carried out in [5] 
and concluded that space vector pulsewidth modulation 
(SVPWM) algorithm gives good performance. The 
generation of switching pulses in SVPWM algorithm is 
given in [6]. Though the SVPWM algorithm gives good 
performance, the complexity involved is more due to 
angle and sector calculations. To reduce the complexity 
involved in SVPWM algorithm, carrier based SVPWM 
algorithm is developed in [7] by adding the offset voltage 
to the reference phase voltages. Hence, to reduce the 
ripples and to obtain the constant switching frequency 
operation of the inverter, SVPWM algorithm is used for 
DTC in the literature [8]-[9]. 

The two-level inverter fed DTC drives are suitable for 
low power applications only. In order to meet the 
medium and high power applications, nowadays 
multilevel inverters are becoming popular. A diode 
clamped three level inverter fed induction motor drive is 
presented in [10]. A detailed survey on various types of 
multilevel inverters is given in [11]. But, as the number 
of levels increases, the complexity involved in the PWM 
algorithms also increases. To reduce the complexity 
involved in the SVPWM algorithm, a simplified 
SVPWM algorithm is presented for three-level inverter in 
[ 12 ]. 

To reduce the complexity, a novel voltage modulation 
algorithm is presented in [13] by using the concept of 
effective time. By extending the same concept, various 
PWM algorithms have been generated for two-level 
inverters in [14]. The concept of the PWM algorithms 
which is presented in [13]- [14] is extended for multilevel 
inverters in [15]. 
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This paper presents a simple generalized scalar PWM 
algorithm for 2-, 3-, 5- and 7-level inverter. The proposed 
algorithm has been developed by using the concept of 
imaginary switching times, which are proportional to the 
instantaneous sampled phase voltages only. Moreover, 
the proposed algorithm does not require the calculation of 
angle and sector information and hence reduces the 
complexity involved in the PWM algorithm. 



T 2 =^-(-V d )*T s =^(-j3*V d ) 

v dc v dc 
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S V bn - -r~V cn = T bn - T a 



2. Simplified SVPWM Algorithm for Two- 
level inverter 

To reduce the complexity involved in the conventional 
SVPWM algorithm and execution time due to longer 
calculations, a simplified SVPWM algorithm has been 
developed using the concept of imaginary switching 
times. In this approach, the actual switching times for 
each inverter leg are deduced directly as a simple form. 
The proposed approach is based on the instantaneous 
values of the reference voltages of a, b and c phases only. 
This method does not depend on the magnitude of the 
reference voltage space vector and its relative angle with 
respect to the reference axis. The transformation from 
two-phase voltages to three-phase voltages can be 
obtained from the stationary frame reference voltages as 
given in (1) 




If the reference voltage vector lies in the first sector as 
shown in Figure 1 , then the actual switching times can be 
deduced as 
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But from Figure 1 it can be observed that 
V q - V r ef cos a an< 3 Vj = - V ref s i n a • Hence, the active 

vector times can be simplified as 
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Figure. 1 Voltage space vectors produced by an inverter 



From (4) - (5), the imaginary switching time periods 
proportional to the instantaneous values of the reference 
phase voltages are defined as 



T = 

± an 



\V dc J 



1 bn 



yVdc J 



v bn 



( 6 ) 



T = 

± cn 



V Vdc J 



Hence, from (4) - (5), the active vector switching times 
Ti and T 2 , if the reference voltage vector falls in sector- 1 
may be expressed as 

(V) 



T _ T T 

1 \ 1 an 1 bn 



T 1 ~ T bn ~ T cn 



Thus, the active voltage vector switching times can be 
represented by the time difference between the imaginary 
switching time periods. In the SVPWM algorithm, when 
the reference voltage vector falls in the first sector, the 
imaginary switching time which is proportional to the a- 
phase ( T an ) has a maximum value, the imaginary 



switching time which is proportional to the c-phase ( T cn ) 
has a minimum value and the imaginary switching time 
which is proportional to the b-phase ( T bn ) is neither 
minimum nor maximum switching time. Thus, in general 
to calculate the active vector switching times, the 
maximum and minimum values of imaginary switching 
times are calculated in every sampling time as given in (8) 

-(9). 



and 



T max= M ax(T an ,T bn ,T cn ) 

T mm = Min{T an ,T bn ,T cn ) 



( 8 ) 

(9) 
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The effective time T e g can be defined as the time 
difference between T max and T min and can be given as 
in (10). 

Teff = ^max “ ^min (^) 



The effective time means the duration in which the 
effective voltage is supplied to the machine terminals. In 
the actual switching instants, there is one degree of 
freedom that the effective time can be located anywhere 
within one sampling interval. To generate actual 
switching pattern which preserves the effective time, the 
zero sequence time is subjoined to the phase voltage time. 
In order to locate the effective time in centre of the 
sampling interval, the zero sequence voltage has to be 
symmetrically distributed at the beginning and end of one 
sampling period. Therefore, the actual switching times 
for each inverter leg can be simply obtained by the time 
shifting operation as below. 

T ga = T as + T offset 

T gb= T bs+ T offset (H) 

T gc = T cs + 7 offset 

To distribute zero voltage symmetrically during one 
sampling period, the offset time T 0 ff set is achieved using 

a simple sorting algorithm, the zero voltage vector time 
duration can be calculated as given in (12). 

T zero = T s ~ T eff ( 12 ) 

And, T m j n +T 0 jf set > T zero / 2 (13) 

Therefore, T offset = T zero / 2 - T mm (14) 

In order to generate symmetrical switching pulse pattern 
within two sampling intervals, when the switching 
sequence is ‘ON’ sequence, the actual switching time 
will be replaced by the subtraction value with the 
sampling time as follows: 

Tga,gb,gc ~T S ~Tga,gb,gc (15) 

As described above, the effective time implies the 
applied time of a certain active vector. Therefore, with 
the effective vector concept, the actual switching time 
can be obtained directly from the stationary frame 
reference voltage without sector identification, effective 
time calculation and recombination. 

To implement the novel space vector PWM, only a 3- 
element simple sorting algorithm is needed. Therefore, 
the calculation efforts of the proposed PWM method is 
greatly reduced as compared with the conventional space 
vector PWM method. 

3. Topology of n-level Diode-Clamped 
Multilevel Inverter 

An inverter which is capable of three-level or more is 
termed as Multilevel Inverter. The topology of n-level 
diode clamped multilevel inverter is as shown in Figure. 
2. In this circuit, the dc-bus voltage will be split into 
suitable number of levels according to the number of 



capacitors. Although each active switching device is only 



required to block a voltage level of 



Vdc 

(n- 1) 



the clamping 



diodes must have different voltage ratings for reverse 
voltage blocking. Assuming that each blocking diode 
voltage rating is the same as the active device voltage 
rating, the number of diodes required for each phase will 
be (n - 1) x (n - 2) . This number represents a quadratic 



increase in n. when n is sufficiently high, the number of 
diodes required will make the system impractical to 
implement. If the inverter runs under PWM, the diode 
reverse recovery of these clamping diodes becomes the 
major design challenge in high-voltage high-power 
applications. 



4. Proposed Generalized PWM Algorithm 
for an n-Level Inverter 

For a 2-level inverter only one offset time is sufficient to 
evaluate the switching instants. Because, in a two-level 
SVPWM algorithm, in a sampling time interval, the 
imaginary switching time which has lowest magnitude 
( ^min ) crosses the triangular carrier first, and causes the 
first transition in the inverter switching state. While the 
imaginary switching time, which has the maximum 
magnitude ( r max ), crosses the carrier last and causes the 
last switching transition in the inverter switching states. 
Thus the switching periods of the active vectors can be 
determined from the T max and Tj^ in a two-level 



inverter scheme. 

But, the SVPWM algorithm for multilevel inverters, 
involves comparing the modulating signals with a 
number of symmetrical level-shifted triangular carrier 
waves for PWM generation. It has been shown that for an 
n-level inverter, n-1 level-shifted triangular carrier waves 
are required for comparison with the modulating 
waveforms. Because of the level- shifted multicarriers, the 
first crossing (termed the first-cross) of the modulating 
wave cannot always be the min-phase. Similarly, the last 
crossing (termed the third-cross) of the modulating wave 
cannot always be the max-phase. Thus, a single offset 
time is not sufficient to centre the middle inverter 
switching vectors, in a multilevel PWM scheme during a 
sampling period Ts. 

The proposed algorithm determines the offset time for the 
PWM generation in the linear modulation region by using 
the concept of imaginary switching times. The proposed 
generalized PWM algorithm presents a simple way to 
determine the time instants at which the three reference 
phases cross the triangular carriers. These time instants 
are sorted to find the offset voltage to be added to the 
reference phase voltages for SVPWM generation for 
multilevel inverters for the entire linear modulation range, 
so that the middle inverter switching vectors are centered 
(during a sampling interval), as in the case of the 
conventional two-level SVPWM scheme. 

For an n-level inverter, the instantaneous imaginary 
switching times are calculated as follows: 



T as = V an x 



VjcKn- 1 ) 



(16) 
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where n represents number of levels. 





The time equivalent of a common mode voltage is termed 
as T offset! which can be defined as 

T offset! ~ — (-^max ^min ) ^ (19) 

where T max and r min are the maximum and minimum 
oiT as’ T bs and T cs . 






Figure. 2 Topology of n-level multilevel inverter 



Then time equivalents of the 
magnitudes can be defined as 


modified phase voltage 


so 

+ 

>5 

II 


(20) 


^hs Tfe T offset 1 


(21) 


Tcs = T cs + T 0 ffset\ 


(22) 



The triangular carriers and the modulating waves, for an 
n-level PWM scheme are shown in Figure. 3a, for n is 
odd, and in Figure. 3b, for n is even. The (n-1) triangular 
carriers are compared with reference phase voltages as 
shown in Figures. 3a and 3b. A carrier index, I, is defined 
to designate the carrier regions in which the reference 
phase voltages lie during the sampling interval under 
consideration. During a sampling interval, the carrier 
indices, I a , I b and I c (which can be from 1 to n-1), for a, b 
and c phases, respectively, are determined depending on 
the carrier region in which the respective phase voltage 
lies. Then the time durations, T a _ cross , T b _ cross and 



^c- cross when n is odd, can be determined for a n-level 



inverter as: 

Ta _ cross = T*s + «/ fl - (n - 1) / 2) * r, (23) 

T b _ cross = + ((h - (« - 1) 1 2 ) * T s (24) 

T c _cross = T cs +(dc ~(n-l)/2)*T s (25) 

When n is even, the time durations, T a _ cross , T b _ cross 
and T c _ cross , can be determined for a n-level inverter as: 

T a _cross =(T s /2)+T* s+ ((I a -(n/2))*T s ) (26) 

T b _ cross ={T s /2)+TZ s +((I b -(n/2))*T s ) (27) 

T c _cross=(T s /2) + T* s +((I c ~(n/2))*T s ) (28) 



in the proposed generalized PWM algorithm, the 
^a— cross > ^b— cross T c — cross are used to centre the 
middle switching vectors as in the case of two-level 
inverters in each sampling time interval. The time- 
duration, at which the reference phases cross the 
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triangular carriers for the first time, second time and third 
time are defined as T first _ cross » ^second _ cross 
T t hi r d cross respectively, in a sampling time interval T s . 
these can be calculated from (29). 



1 first _cross 

T 

L sec ond _ cross 



- min' 



Across) 

— mid{r x cross ) 



(29) 



' third _ cross 



= max 



^fx_cross\ % a,b,i 





(b) 

Figure. 3 Triangular carrier and reference signals (a) n-level PWM 
scheme where n is even (b) n-level PWM scheme where n is even 



During one sampling interval Ts, the time durations 
T first _ cross * ^sec ond _ cross ^ third _ cross directly 

decides the switching times for different inverter voltage 
vectors. An offset time, T offset2 is added to T first cross , 

Lee ond _ cross and T third _cross to centre the middle 
vector. The offset time, T 0 ^ set2 is determined as follows: 

The time duration for the middle inverter switching 
vectors, T middle , is given as 

middle ~~ ^ third _ cross ~T first _ cross (30) 

The time duration for start and end vector is 

T 0 = T s ~ T middle d) 

The time duration for start vector is given by 

T 0 /2 = T first _ cross + T 0 ff se t2 ( 32 ) 

Therefore, T oJfset2 =T 0 /2-T first _ cross (33) 



By adding the time, T ( ffet2 to Tf r st _cross > 
■^see ond _ cross and T third _ cross gives the inverter leg 
switching times T ga , T gb and T gc for phases a, b and c, 
respectively as (34)-(36). 

T ga =T a_ cross + T offset! ( 34 ) 

Tgb ~ T b _ cross + T 0 jf set 2 (35) 

Tgc -T c _ cross +T 0 jfset! (36) 

5. Results and Discussion 

To validate the proposed generalized scalar PWM 
algorithm, several numerical simulation studies have 
been carried out and results are presented. The details of 
the induction motor, which is used for simulation studies 
are as follows: 

A 3-phase, 4 pole, 4kW, 1200rpm induction motor with 
parameters of R s = 1.57D, R r = 1.21^, L s = L r = 0.17H, 
L m = 0.1 65H and J = 0.089Kg.m 2 . 

The steady state simulation results for 2-, 3-, 5-, and 7- 
level inverter fed DTC-IM drive are shown from Figure. 
4 to Figure. 19. From the simulation results it can be 
observed that as the number of levels increases the 
harmonic distortion will be decreased. Also, it can be 
concluded that the proposed PWM algorithm is efficient 
with reduced complexity. The total harmonic distortion 
(THD) values of line currents and voltages are given in 
Table- 1. 



TABLE I 

COMPARISON OF THD IN LINE CURRENT FOR VARIOUS LEVELS 



Level 


% THD in line 
current 


% THD in line voltage 


2 level 


9.84 


75.96 


3 level 


4.06 


38.10 


5 level 


2.23 


17.25 


7 level 


1.37 


12.17 




Time (-sec) 

Figure. 4 Simulation results of 2-level SVPWM based DTC: steady- 
state plots of speed, torque, currents and flux at 1200 rpm 
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0.35 0.4 0.45 

Time (sec) 



Figure. 5 Simulation results of 2-level SVPWM based DTC: steady 
state plots of phase and line voltages 




Figure. 6 Harmonic Spectrum of line current for 2-level SVPWM 
based DTC along with THD 




Figure. 7 Harmonic Spectrum of line voltage for 2-level SVPWM based 
DTC along with THD 
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Figure. 8 Simulation results of 3-level SVPWM based DTC: steady-state 
plots of speed, torque, currents and flux at 1200 rpm 
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Figure.9 Simulation results of 3-level SVPWM based DTC: steady 
state plots of phase and line voltages 




Figure. 10 Harmonic Spectrum of line current for 3-level SVPWM 
based DTC along with THD. 




Figure. 1 1 Harmonic Spectrum of line voltage for 3-level SVPWM 
based DTC along with THD. 
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Figure. 12 Simulation results of 5 -level SVPWM based DTC: steady- 
state plots of speed, torque, currents and flux at 1200 rpm 




Figure. 13 Simulation results of 5 -level SVPWM based DTC: steady 
state plots of phase and line voltages 




Figure. 14 Harmonic Spectrum of line current for 5-level SVPWM 
based DTC along with THD 
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Figure. 15 Harmonic Spectrum of line voltage for 5-level SVPWM 
based DTC along with THD 
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Figure. 16 Simulation results of 7-level SVPWM based DTC: steady- 
state plots of speed, torque, currents and flux at 1200 rpm 




Figure. 17 Simulation results of 7-level SVPWM based DTC: steady 
state plots of phase and line voltages 




Figure. 18 Harmonic Spectrum of line current for 7-level SVPWM 
based DTC along with THD 




Figure. 19 Harmonic Spectrum of line voltage for 7-level SVPWM 
based DTC along with THD 



6. Conclusion 

A simple and novel generalized SVPWM algorithm is 
presented for an n-level inverter fed direct torque 
controlled induction motor drive. The proposed algorithm 
has been developed by using the concept of imaginary 
switching times, which are proportional to the sampled 
reference phase voltages only. Hence, the complexity 
involved id less when compared with the classical PWM 
algorithms. To validate the proposed algorithm, the 
simulation studies have been carried out and results are 
presented. From the simulation results, it can be 
concluded that the proposed algorithm gives good 
performance with reduced complexity. Moreover, as the 
number of levels increases, the harmonic distortion of 
line current and voltages also decreases. 
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Abstract 

Ce travail synthetise revolution de legalisation 
classique vers legalisation par reduction de canal. II 
s’agit d’utiliser les statistique derdre superieure a 
deux a savoir le calcul du maximum de kurtosis pour 
retrouver le filtre TEQ (Time-domain EQualizer) op- 
timal. Le delai de synchronisation dans l’algorithme 
du MSSNR (Maximum Shortening Signal to Noise 
Ratio) est retrouve par la meme technique. Une 
comparaison entre le filtre optimal retrouve par 
d’autres methodes et celui determine par la methode 
la notre. La fonction statistique gamma a reduire de 
maniere consequente le delai de convergence apres 
son introduction dans l’algorithme. 1 

Keywords: TEQ (Time-domaine EQalizer), 

High Order Statistics, kurtosis, DMT, fonction 
Gamma. 

1 Introduction 



optimales et/ou adaptatives [12] realisent le meme 
travail avec moins de performances. La plupart 
des approches par modele TEQ (Time domaine 
Equalizer) ne sont pas adaptatives, elles necessitent 
l’apprentissage ou une estimation du canal avec des 
calculs assez lourds a mettre en oeuvre. L’algorithme 
ameliore est celui de MERRY (Multicarrier Equali- 
sation by Restoration of RedundancY) [6, 11]. Nous 
proposons une methode basee sur les HOS, dite par 
maximum du kurtosis, et nous calculons egalement le 
delai optimum ainsi que le filtre optimal lui correspon- 
dant. 

Cet article est organise de la fagon suivante: Section 
2 met en evidence l’algorithme d’egalisation multi- 
porteuse par retablissement de redondance.La fonc- 
tion cout sera donnee en Section 3. Notre contribu- 
tion sera mentionnee en Section 4. Des simulations 
et discutions des resultats fairront l’objet de la partie 
theorique en Section 5. Les resultats importants ainsi 
que les futurs travaux envisages seront donnes en con- 
clusion. 



L ’evolution de legalisation classique vers legalisation 
par reduction de canal a fait l’objet de nombreux 
travaux dans la literature. Nous nous sommes in- 
teresses a l’introduction des HOS dans ce damaine. 
Nous appliquons le calcul du maximum de kurtosis 
pour retrouver le filtre TEQ (Time-domain EQual- 
izer) optimal ainsi que le delai de synchronisation dans 
l’algorithme du MSSNR (Maximum Shortening Sig- 
nal to Noise Ratio). Nous comparons le filtre optimal 
retrouve par d’autres methodes et celui determine par 
la methode du maximum du kurtosis. Nous arrivons 
par utilisation des fonctions statistiques notamment 
gamma, a reduire de maniere consequente le delai de 
convergence au prix de la stability de 1’ algor it hme. 
Des methodes retrouvent le canal et calculent 
l’egaliseur optimal [8]. D’autres sont dites sous- 



1 Ce travail a ete realise au departement des TCSN, FSNV. 



2 Algorithme d’egalisation mul- 
tiporteuse par retablissement 
de redondance 

L’algorithme MERRY (Multicarrier Equalization by 
Restoration of RedundancY) exploite la redondance 
introduite dans les donnees a travers le prefixe cy- 
clique, voir figure (2) [6]. C’est une approche basee 
sur la minimisation de l’esperance mathematique du 
carre de la difference cyclique qui est la difference en- 
tre des valeurs de separees par un bloc de donnees. 



y( 2) = h()s{2) + his(l) + h 2 s( 0 ) + h^s(-l) + ft-4s(— 2) 

= ho<s(10) + h\s{9) + [^ 2 ^( 0 ) + h%s( 1) -h h^s^— 2)] 
7/(10) = ho<s(10) + his(9) + [^2^(8) + h^s{7) + ^45(6)] (1) 
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avec 5 G {0, • • • , M — 1}, S etant le delais desire (c’est 
le coefficient de synchronisation) . Le gradient stochas- 
tique de (5) donne l’algorithme de reduction de canal 
adaptatif aveugle MERRY (Multicarrier Equalization 
by Restoration of RedundancY) [MAR 05], avec: Pour 
un S donne, le symbole k = 0, 1, 2, • • • , 



Figure 1: Rajout du CP a la trame de donnees 



La figure 2 fait appel a l’equation 1. Ainsi, 
les N + P donnees transmises deviennent peri- 
odiques. Cette periodicite est necessaire pour main- 
tenir l’orthogonalite des multiporteuses [10]. Les 
derniers symboles sont identiques aux premiers et 
nous pouvons ecrire : 

s(Mk + i) = s(Mk + i + N ), ,N} (2) 

ou M = P + N est la duree de tous les symboles et 
k est l’indice de chaque symbole. La figure 2 montre 
un exemple ou Ton prend TV = 8, P = 2 et M = 
10 correspondant a k = 0. Les donnees regues sont 
obtenues a partir de selon la relation suivante : 

L c 

r(Mk + i) = ^ ci • s(Mk + i — l) + n(Mk + i) (3) 
z=o 

Les donnees egalisees sont obtenues a partir de suivant 
cette equation : 



Lf 

y(Mk + i) = fj ■ r(Mk + i— j) (4) 
j=o 

La relation donnee par (1) est deterioree a cause du 
canal lors de la reception des donnees, car les ISI qui 
affect ent le CP sont different s de ceux qui affect ent 
les P derniers elements du symbole. Considerons 
F exemple de la figure , Les elements transmis 2 et 10 
sont identiques. Cependant, au niveau du recepteur, 
F interference des elements avant F element 2 ne sont 
pas toutes egales a celles qui sont avant Felement 10. 
Si on considere /i 2 , hs et Annuls, alors y( 2) = y{ 10) car 
5(2) = s(10) et s( 1) = 5(9). II est clair que pour avoir 
il faut forcer h 2 , ^3 et h 4 a zero (elements se trouvant 
entre crochets dans Fequation 1), cette egalite etant 
prise dans le sens des moindres carres. Nous rendons 
de cette fagon le canal effectif de taille equivalente a 
celle du CP. 



3 Algor it hme MERRY et sa 
fonction cout 



La fonction cout de l’algorithme MERRY est : 



J = E 



| y(Mk + P + 5) 



y (Mk + P + +N + 5) | 2 

( 5 ) 



r k) = r (Mk + P + <5) — r (Mk + P + N + 5) 

e{k) = f T (fc)rfc) 
i(k + 1) = i T (k) — /^e(k)ik) 

f (k + 1) 



f (k + 1) = 



i(k + 1) 



(6) 



avec r(i) = [r(i), r(i — 1), • • • ,r(i — Lf] T , et ( * ) est 
le conjugue complexe. Pour eviter la solution triviale 
(f = 0 ) nous devons prendre ||f|| = 1. L’algorithme 
represente par (6) donne le vecteur propre minimal de 
la matrice A = E [f , r H ] . Soit C la matrice canal de 
convolution autrement dit h = Cf et soit C wa nobtenue 
par suppression de S lignes des S + P — 1 de C. 
Si l’entree s(k) est blanche, alors A = C^ a ^,C wa n- 
L’energie a l’exterieur de la fenetre du canal effectif 
plus le gain du bruit est donnee comme suit : 



J s = 2ai E l^'| 2 ]+ 2fTR « r ( ? ) 

\J=0 j=5+p 



Dans le cas ou la reduction du canal utilisee est au- 
todidacte et non adaptative, alors la matrice A = 
E [f , r H ] peut etre estimee a partir des donnees. Le 
vecteur propre correspondant a la valeur propre min- 
imale peut etre calcule, ce qui fournira une technique 
d’initialisation pour eviter les modes de convergence 
lents. 

Reduction optimale du canal 

La reduction de canal optimale utilise un algorithme 
base sur le calcul des valeurs propres et des vecteurs 
propres pour generer les coefficients du filtre reduit 
optimal. La technique consiste a determiner le filtre 
TEQ qui maximise le rapport de reduction du signal 
sur le bruit dans une fenetre. Pour ce faire, il faut 
minimiser l’energie de la reponse impulsionnelle du 
canal effectif, a l’exterieur de la fenetre contenant les 
(P+1 ) echantillons. Nous reprenons C^ n et C wa u tels 
que definies dans [8], alors h w i n = C w i n f represente le 
canal effectif situe dans la fenetre de taille ( P + 1). 
Cependant h wa n = C wa nf donne le canal effectif en 
dehors de cette fenetre. Le probleme de conception 
du TEQ a base du MSSNR peut etre decrit comme 
suit [8, 1, 2, 4, 9]: 



min (f T Af) sous la contrainte f T Bf = 1 (8) 



ou A et B sont des matrices reelles [MAR 02] , 
symetriques Lf x Lf avec : 



A — C wa n , Cwalli B — G w in , 



( 9 ) 
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Elies represented, respectivement, les energies corre- 
spondant a C wa uet C win ou : 



hlallKall={ T CZ all C{={ T M ( 10 ) 



uT -r fTpT p r fTp>r 

LL win Li win — 1 ^ win-win 1 ~ 1 01 



(ii) 



La resolution de l’equation (8) mene au TEQ qui 
donne la solution du probleme de vecteurs propres 
generalise [10], nous avons par consequent : 



Bf = AAf 



( 12 ) 



ou A est la valeur propre generalisee correspondant 
au vecteur propre generalise f, qui est la solution de 
l’equation (12). 

C = (v 7 ^) _1 A (\/s) _1 (13) 



L’ inverse de la mat rice A est retrouve par factorisation 
de Cholesky. Les TEQ sont alors donnes par : 




ou q est le vecteur propre correspondant a la valeur 
propre A. II est montre aussi [8] que l’algorithme 
converge globalement si q = q m i n est le vecteur pro- 
pre minimal correspondant a la valeur propre mini- 
male A min • C’est ainsi que le Litre TEQ optimal 
foptest retrouve [3]. L’expression du SNR reduit opti- 
mal (Shortening SNR ou SSNR) est donnee par : 

SSNR opt = 10 log ( = 10 log ( -?—) 

\ ^opt^opt J \ Amin J 

SSNR op t = —10 log (Xmin) (15) 

En somme, Lalgorithme optimal utilise la matrice C, 
dont les elements sont calcules a partir de H wa n et 
H w i n . Le SSNR optimal est retrouve par calcul direct 
des valeurs propres minimales de la matrice C. Les 
coefficients du TEQ optimal sont finalement calcules 
par transformation lineaire du vecteur propre unitaire 
associe a la valeur propre minimale. 



4 Modification de l’algorithme 

4.1 Le filtre optimal 

Cette procedure debute a partir de tous les filtres re- 
duits f n ayant ete trouves a partir de l’equation (14), 
utilisant l’algorithme du MSSNR [8]. Nous proposons 
de choisir le filtre SCE (Shortened Channel Equalizer) 
optimal, non pas en comparant les energies de tous les 
filtres trouves tels que abordes par Martin et al. dans 
[7], mais plutot en calculant le kurtosis de chaque fil- 
tre. Nous avons ainsi : 



Reduire le canal revient a trouver un filtre tel que sa 
sequence de sortie ait la meme variance et le meme 
kurtosis maximum en valeur absolue qu’a l’entree du 
filtre. Le filtre optimal sera celui qui maximisera la 
quantite donnee par l’equation (19). Le delai de syn- 
chronisation sera la valeur de l’indice qui correspondra 
a ce filtre, tel que: 

k = arg{max(g n )} (18) 

n 

4.2 Calcul du pas de convergence par 
kurtosis 

Le kurtosis normalise d’une variable aleatoire s ayant 
comme valeur moyenne s est defini par le rapport du 
moment d’ordre quatre 774 au carre du moment d’ordre 
deux 77 I f° us deux centres. Les HOS (Higher Order 
Statistics) definissent le kurtosis comme suit [LAC 97]: 

«4(S) = ^ (19) 

V2 

ou, si est E Toperateur esperance mathematique: 

V2 = E {(s - s) 2 } (20) 

t ?4 = £((s-s) 4 j (21) 

L’exces (kurtosis standard) est defini par l’equation 
(19) moins trois. II s’annule lorsqu’il s’agit d’une dis- 
tribution normale [3, 5]. D’ou tout l’interet a intro- 
duire ce coefficient. En effet, sachant que le bruit ad- 
ditif est dans la plupart des cas de type Gaussien, il 
est bien connu que les statistiques d’ordre superieur ne 
sont pas affectees par le bruit. Nous proposons dans 
ce papier une methode iterative pour le calcul de fi. En 
effet, le pas de convergence divise par le kurtosis nor- 
malise donne par l’equation (19) de la sequence regue a 
l’entree de l’egaliseur permet de tenir compte a chaque 
iteration de la sequence a l’entree de l’egaliseur r(i). 
Ainsi, le pas de convergence du nouvel algorithme est 
calcule a partir d’un fiQ divise a chaque iteration par 
le kurtosis de sorte a avoir: 

Mo 

W+1 “k 4 ( r(*)) 

Par ce precede, la mise a jour du pas de convergence 
depend des donnees regues (c’est a dire des conditions 
de transmission du canal). Deplus, nous utilisons 
la fonction de densite de probability Gamma connue 
dans la litterature pour la determination du pas de 
convergence. Les result at s obtenus vont maintenant 
etre discutes dans la partie consacree aux simulations. 

5 Simulations 



y„ = h * f n 

gn = «4(y n ) 



^ 5 ) Les donnees que nous utilisons lors de la simulation 
sont les huit CSA loop (Carrier Serving Area) [7], util- 
(17) isees dans le standard ADSL. Nous utilisons le CSA 1 
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lors de nos simulations. Le CP est de taille, la taille 
de Tegaliseur est egale a 16 [9]. La figure 2 represente 
les reponses impulsionnelles du canal reelle et celle 
reduite par MSSNR a partir de Talgorithme MERRY 
optimal. La figure3 represente les reponses impulsion- 




Figure 2: Canal reel et canal reduit par MSSNR 

nelles du canal reelle et celle reduite par la methode 
du maximum du kurtosis. Afin de pouvoir comparer 
les result at s donnes par les deux methodes, nous nous 
proposons de derouler 1’ algor it lime de MERRY et 
celui que nous proposons. Le pas de convergence etant 




Figure 3: Canal reel et canal reduit par maximum du 
kurtosis et la fonction gamma 



determine par notre technique d’iteration dependant 
de la sequence de donnees ou de fagon classique. Dans 
tous les cas nous choisissons un SNR egal a 40dB. 

Les figures (2 et 3) donnent les reponses impul- 
sionnelles reelle et reduite dans les deux cas (a savoir 
la methode du MSSNR et celle du kurtosis conjuguee 
a la fonction Gamma. Les deux figures sont tres sem- 
blables mis a pars quelques legeres differences. La 
figure 4 donne la fonction cout de MERRY ainsi que 
le debit. Nous remar quons que la convergence est ef- 
fective a partir de la 500eme iteration avec un debit 
fluctuant entre 3 et 4Mb/s. La figure 5 represente la 
fonction cout et le debit faisant appel aux statistiques 
d’ordre superieur HOS et la fonction Gamma. Nous 
remarquons que la convergence est effective a partir 
de la 180eme iteration avec un debit fluctuant entre 3 
et 4Mb/s. 





Figure 4: Courbe de convergence par MERRY optimal 





Figure 5: Courbe de convergence par kurtosis et la 
fonction Gamma 



II est a remarquer que les delais optimaux sont 
differents de ceux retrouves dans la litterature [6]. 



6 Conclusion 

Nous avons ameliore les approches adaptatives de re- 
duction de canal aveugle. En effet, une nouvelle meth- 
ode de calcul du filtre optimal par maximisation du 
kurtosis a ete proposee. Le delai de synchronisation 
a pu etre calcule. Les resultats trouves ont ete com- 
pares avec ceux donnes dans la litterature. A travers 
les HOS (High Order Statistics) conjugue a la fonction 
Gamma le temps de convergence de Lalgorithme a ete 
nettement reduit. II est a signaler que la stability de 
la fonction cout et du debit est meilleure avec la meth- 
ode de MERRY. La stability des resultats fera l’objet 
de travaux futurs. Tous nos travaux de simulation ont 
ete realises sous Matlab. 



Remerciements 



Nous remercions le Pr J. M. Rouvaen du laboratoire 
IEMN de l’universite de Valenciennes en France pour 



50 



ICGST-ACSE Journal, ISSN 1687-481 1 , Volume 12, Issue 2, Delaware, USA, October 2012 



son aide. Nous remer cions aussi, l’universite de Texas 
AT, Austin aux USA pour avoir mis a notre disposi- 
tion les donnees standards CSA. 



References 

[1] N. Al-Dhahir. FIR Channel Shortening Equal- 
izers for MIMO ISI Channels. IEEE Trans, on 
comm , 49(2):213-218, 2001. 

[2] G. Arslan, B. Lu, L. D. Clark, and B. L. 
Evans. Iterative Refinement Methods for Time 
Domain Equalizer Design. EURASIP Processing , 
2006(1) :83-95, 2005. 

[3] I. Chahed, H. Bellahsene, J. M. Rouvaen, and 
M. Djeddi. Selection of the order of ARMA 
model using Kurtosis minimisation. Application 
to blind deconvolution of seismic data. AMSE 
proceding , 47(2):77-89, 2004. 

[4] D. D. Falconer, S. L. Ariyavisitaul, A. Benyamin- 
Seeyar, and B. Eidson. Frequency domain equal- 
ization for single-carrier broadband wireless sys- 
tems. IEEE Comm. Mag., 40(4):58-66, 2002. 

[5] J.L. Lacoume, P.O. Amblard, and P. Comon. 
Statistiques d’Ordre Superieur pour le Traitement 
du Signal. Masson, 1997. 

[6] R. K. Martin, J. Balakrishnan, W. A. Sethares, 
and C. R. Johnson. A Blind Adaptive TEQ for 
Multicarrier Systems. IEEE Signal Processing 
Letters , 9(ll):341-343, 2002. 

[7] R. K. Martin and C. R. Johnson. Adaptive 
Equalization: Transitioning from single-carrier to 
multicarrier systems, in Proc. IEEE Signal mag., 
22 (6) : 108—122, 2005. 

[8] P. J. W. Melsa, R. C. Younce, and C.E. Rohrs. 
Impulse Response Shortening for discrete Mul- 
titone Transcievers. IEEE Trans, on comm., 
44(12): 1662-1672, 1996. 

[9] H. Minn, N. Al-Dhahir, and Y. Li. Opti- 
mal Training Signals for MIMO OFDM Chan- 
nel Estimation in the Presence of Frequency Off- 
set and Phase Noise. IEEE Trans, on comm., 
54(10):1754-1759, 2006. 

[10] T. Pollet, M. Peeters, M. Moonen, and L. Van- 
dendorpe. Equalization for DMT-Based Broad- 
band Modems. IEEE comm, maq ., 38(5): 106- 
113, 2000. 



Biographies 



Dr. Bellahsene Hocine a regu 
son diplome d’ingenieur en communica- 
UpL I tion en 1993 et le diplome de Magistere 
1999 a l’universite de Tizi-ouzou 
(Algerie). II prepara sont diplome de 
doctorat au laboratoire IEMN labora- 
tory (Valenciennes, France) de janvier 
| 2005 a juillet 2007. II regu son doctorat 

en juin 2009. II est actuellement maitre de conference a 
l’universite de Bejaia (Algerie). Son domaine d’interet 
porte sur legalisation aveugle dans les systemes de 
coommunications, les statistiques d’ordre superieur et les 
techniques MIMO pour les communications numeriques. 
mail: bellahsene ho@yahoo.fr. 

Fatani Imade Fahd Eddine est 

actuellement en post doc a l’univeriste 
de Valenciennes en France. II a obtenu 
son diplome de docteur en telecom- 
munication a la meme universite 
en septembre 2010. Actuellement il 
travail a l’ecole nationale des sciences 
appliquees de khouribga au Maroc. 
Son domaine d’interet porte sur les techniques de codage 
conjoint source canal, les techniques MIMO pour les 
transmissions et les transmission video H 264 dans le 
contexte MIMO. Mail: fatani.imade@gmail.com 




[11] A.H. Sayed. Fundamentals of adaptive filterinq. 
John Wiley, 2003. 

[12] B. Widrow, J. M. McCool, M. G. Larimore, and 
C. R. Johnson. Stationary and Nonstationary 
Learning Characteristics of the LMS adaptive fil- 
ter. Proc. IEEE, 64(8):1151-1161, 1976. 



51 



